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1  Introduction 

The  primary  goal  of  this  research  is  to  devdop  a  mesoscale  model  capable  of  predicting 
rc^oos  of  convective  outbreak.  The  modd  is  to  be  initialized  by  a  simple  sounding(s) 
and  during  the  data  assimilation  period  sin^e  Doppler  radar  data  are  to  be  assimilated 
into  the  model  fields.  The  model  should  be  capable  of  producing  forecasts  of  convective 
development  on  a  mini-casqmter  system.  The  modd  should  cover  a  100  km  by  100  km 
domain  with  a  minimum  resolution  ol  5  km.  The  benchmark  mini-computer  system  used 
for  modd  developniant  is  the  DEC  Micro- Vax  n.  The  basic  modd  used  for  this  study  is 
the  Colorado  State  University  Regional  Atmospheric  Modeling  System(RAMS).  The  overall 
featoces  of  RAMS  ore  described  in  Appendix  A.  hi  the  fiillowing  sections  we  describe  the 
various  vecakos  of  RAMS  that  have  been  developed  for  this  project,  the  methods  of  four- 
cBmsttdooal  asdmilation  of  sin^  Doppler  radar  data,  preliminary  testing  and  evaluation 
procedures  and  results.  We  then  describe  the  real  data  cases  used  to  test  the  modd,  the 
methods  of  single  Doppler  wind  retrieval,  and  the  results  hf  testing  the  modd  in  several 
cases.  We  condude  with  an  overall  summary  of  results  and  recommendations  for  future 
research. 
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2  Model  development 


Two  versions  of  the  Regional  Atmoepheric  Modelling  System  (RAMS)  were  employed  for 
this  study.  These  were  (1)  a  simplified  version  of  the  full  hydrostatic  model  with  fewer 
grid  points  than  usually  used  in  a  research  mode  and  (2)  a  mixed  layer  modd  developed 
for  this  project,  similar  to  the  modds  of  Lavoie  (1972)  and  Keyscr  and  Anthes  (1977). 
The  two  modd  versioos  were  developed  in  paralld  as  it  was  recognised  from  the  beginning 

the  project  that  each  configuration  had  advantages  and  disadvantages  adien  compared 
to  each  other.  If  execution  speed  were  not  a  factor,  the  full  modd  would  obvioudy  have 
been  the  choice.  One  of  the  advantages  of  the  full  modd  is  it  incorporates  far  fewv 
approximations  than  the  mixed  layer  naodeL  However,  the  full  modd  be 

enou^  to  attain  one  ci  the  stated  goals  of  this  study  which  is  to  produce  a  forecast  in 
1/6  real  time  on  a  MicroVAX  n  class  con^mter  system,which  is  rated  at  about  1  million 
instructions  per  second  (MIPS),  and  to  achieve  this  speed  with  the  horisontal  domain  setup 
of  100  km  by  100  km  total  sise  and  a  6  km  resolution.  The  mixed  layer  modd,  on  the  other 
hand,  has  the  advantages  of  computer  speed  and  dmplicity.  These  modd  versions  will  be 
described  below. 

A.  The  Hydrostatic  FVill  Numerical  Prediction  Modd 

The  fiill  modd  version  uses  the  cooq>lete  hydrostatic  modd  equation  set  from  RAMS. 
This  set  indudes  prognostic  equations  for  the  u  and  v  wind  components,  potential  temperw 
ture,  and  vapor  mixing  ratio  (as  a  pasdve  tracer  only)  and  diagnostic  equations  for  vertical 
vdodty  and  pressure.  The  simplifications  necessary  for  this  project  indude  the  choice  of 
the  simplest  numerical  schemes  and  parameterisations  that  would  perform  adequatdy.  In 
addition,  a  couple  of  the  schemes  needed  to  be  modified  for  the  limited  height  of  the  domain. 
Some  of  these  simplifications  and  modifications  are: 

1.  The  number  of  vertical  levcb  needed  to  be  reduced.  This  limited  the  vertical  extent 
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of  the  atmosphcTe  which  could  be  modelled.  A  domain  of  about  4  km  depth  with  6 
layers  was  first  used.  The  vertical  grid  spacing  was  400  m  near  the  ground,  stretched 
to  750  m  near  the  model  top.  Also,  the  number  of  levels  in  the  soil  model  was  reduced 
from  the  usual  11  to  4. 

2.  The  simplest  torms  of  the  model  numerical  schemes  needed  to  be  used.  These  included 
the  second-order  advection  scheme  rather  than  a  sixth  order,  a  constant  pressure  top 
boundary  condition  for  the  hydrostatic  equation,  and  a  radiative  lateral  boimdary 
condition  with  a  specified  phase  velocity. 

3.  The  diflusion  calculations  were  modified  to  reduce  the  munber  of  calcxilatlons  involved. 
For  instance,  only  the  vertical  shear  was  used  to  compute  the  deformation  rather  than 
the  three  dimensional  shear. 

4.  The  radiative  parameterizations  needed  to  be  modified  to  account  for  the  atmosphere 
above  the  model  domain.  In  the  resetfch  model,  the  domain  usually  extends  well  into 
the  stratosphere  which  will  account  for  almost  all  of  the  mass  of  the  atmosphere.  Since 
the  domain  in  this  case  covers  only  the  lower  portion  of  the  atmosphere,  especially  the 
long  wave  radiation  needed  to  have  information  fipom  the  remainder  of  the  atmosphere. 
The  long  wave  parameterization  was  modified  to  use  the  thermal  information  &om 
the  input  sounding  along  with  the  actual  model  domain. 

6.  Tba  convective  parameterization  also  needs  information  from  all  of  the  troposphere 
and  the  lower  stratosphere.  This  scheme  was  also  modified  to  use  information  from 
the  input  sounding. 

This  model  configuration  was  run  successfully  In  several  situations.  The  convective 
boundary  layer  was  handled  well  and  realistic  slope  flows  were  developed  with  a  ridge-plain 
topography.  The  goal  of  1/6  real  time  was  not  quite  realized;  the  simplified  model  did 
achieve  1/5  real  time  on  a  MicroVAX  II. 
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However,  problexiu  began  to  arise  when  the  model  boundary  layer  velocity  fields  were 
nudged  to  several  generic  ‘‘observed”  fields  meant  to  simulate  Doppler  observations.  As  an 
exaunple  of  these  problems,  consider  the  case  of  a  low  level  stationary  convergence  line  th2it 
is  homogeneovis  in  the  north-south  direction.  In  the  east-west  direction,  the  wind  changed 
from  5  m/s  to  -5  m/s  over  a  20  km  region  in  the  center  of  the  100  km  domain. 

Unlike  the  mixed  layer  model,  the  full  model  has  to  satisfy  the  various  physical  con¬ 
straints  imposed  by  the  equations,  such  as  mass  continuity.  When  a  convergent  field  is 
imposed  at  low  levels,  the  model  response  is  to  increase  pressure  in  the  center  of  the  do¬ 
main  and  force  an  upward  vertical  motion.  The  high  pressure  forces  a  compensating  outflow 
above  the  inflow  layer  with  subsidence  occurring  on  a  larger  scale  than  the  convergent  re¬ 
gion.  The  full  model  has  problems  with  numerical  stability  and  mass  trending  if  dl  the 
parts  of  the  circuUtion  (inflow,  outflow,  subsidence)  are  not  included  in  the  model  domain. 
In  the  model  configuration  required  by  the  computer  speed,  this  was  the  case.  The  domain 
was  not  deep  enough  to  adequately  resolve  the  outflow  and  not  wide  enough  for  the  sub¬ 
sidence  to  be  contained  within  the  domain.  The  mean  domsun  pressrire  increased  to  the 
point  where  numerical  instability  resulted. 

Several  changes  could  be  made  to  the  model  configuration  for  the  nudging  to  maintain 
stability.  More  grid  points  can  be  added  both  horizontally  and  vertically  to  increase  the  size 
of  the  domain.  Obviously,  the  1/6  real  time  goal  on  a  Micro  VAX  II  would  not  be  realized 
as,  along  with  the  increased  number  of  grid  points,  the  timestep  length  would  also  need  to 
be  lowered  since  the  phase  velocity  of  the  modelled  external  gravity  wave  is  dependent  on 
the  height  of  the  model  top.  Another  option  is  to  lower  the  resolution  so  as  to  expand  the 
area  of  the  domain,  but  with  the  expected  loss  of  detail  of  the  modelled  circrilations.  In 
the  validation  simulations  described  below,  the  5  km  resolution  is  retained,  but  more  grid 
points  were  added.  For  the  resd-data  simulations,  the  model  was  used  on  a  100  km  by  100 
km  domain  to  compare  it  with  the  mixed-layer  model.  The  domain  for  both  the  real-data 
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and  validation  simulations  was  increased  to  about  6  km  vertically  by  adding  2  additional 
vertical  levels.  This  model  configtiration  executed  in  about  I/4  retd  time. 

It  appears  that  to  attain  1/6  real  time  forecasts  on  a  Micro  VAX  II  speed  machine  with 
a  100  by  100  km  domain  and  a  5  km  resolution,  the  mixed  layer  version  of  RAMS  must  be 
used.  However,  if  a  faster  machine  is  available,  a  more  complete  model  version  ctin  be  used 
which  has  many  fewer  approximations  than  the  mixed  layer  model. 

B.  The  \Gxed-Layer  Model 

The  mixed-layer-modcl  (MLM)  version  of  RAMS  was  developed  from  similar  models 
described  in  Lavoie  (1972),  Keyser  and  Anthes  (1977),  and  Anthes  et  td.  (1982).  The  latter 
two  references  have  been  followed  most  closely  in  designing  the  RAMS  version.  A  prognostic 
equation  for  moisttire,  which  was  only  found  in  the  first  reference,  has  been  included  here 
also. 

The  purpose  of  a  MLM  is  to  simplify  the  full  3- dimensional  model  equation  set  by 
taking  into  accoimt  certain  gross  features  in  the  vertical  structure  of  the  atmosphere.  Ob¬ 
viously,  this  limits  the  applicability  of  the  model  to  those  problems  having  approximately 
the  vertical  structxire  assmned  in  this  model.  The  assumed  structure  in  particular  consists 
of  a  superadiabatic  surface  layer  50  or  so  meters  deep,  a  neutral  mixed  layer  above  having 
f/Yntt^MTit  vertical  profiles  of  momentum,  potential  temperature,  and  water  mixing  ratio,  a 
possible  inversion  of  zero  vertical  thickness  capping  the  mixed  layer,  and  a  stable  layer  above 
them  aU.  Some  level  within  the  upper  stable  layer  is  called  the  “undistmbed  level”  and  is 
assumed  to  be  the  highest  level  (as  a  function  of  time)  to  which  any  perturbations  generated 
in  the  layers  below  have  penetrated.  The  best  known  atmospheric  situation  which  is  well 
described  by  this  structure  is  the  developing  convective  planetary  boundary  layer. 

Under  the  assumption  of  this  vertical  structure,  six  prognostic  variables  fully  describe 
the  time  varying  atmospheric  state  at  all  heights  below  the  undisturbed  level.  They  are  (1) 
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the  potential  temperature  of  the  mixed  layer,  (2)  the  total  water  mixing  ratio  of  the  mixed 
layer,  (3)  the  westerly  wind  component  in  the  mixed  layer,  (4)  the  southerly  wind  component 
in  the  mixed  layer,  (5)  the  potential  temperature  immediately  above  the  inversion  capping 
the  mixed  layer,  and  (6)  the  height  of  the  mixed  layer  top.  The  field  values  applying  at  the 
undisturbed  level  H  are  obtained  from  the  initial  soun  ling,  and  all  quantities  are  sissumed 
to  vary  linearly  between  H  and  the  mixed-layer  top.  H  itself  is  diagnosed  from  the  mixed- 
layer  height,  and  vertical  velocities  are  diagnosed  from  the  continuity  equation  based  on 
the  mixed-layer  horizontal  wind  components,  the  base  state  density  field,  and  the  s\irface 
topography. 

The  assumptions  about  the  vertical  structure  of  the  atmosphere  permit  the  full  3- 
dimensional  primitive  equations  to  be  integrated  vertically  and  thereby  reduced  to  two 
dimensions.  Thus,  only  a  single  vertical  level  needs  to  be  predicted.  The  prognostic  equa¬ 
tions  for  the  mixed-layer  variables  include  the  forcing  terms  for  horizontal  advection,  hor¬ 
izontal  diffusion,  the  horizontal  pressure  gradient  force,  the  Coriolis  force,  surface  fluxes, 
entrainment  fluxes  at  the  mixed  layer  top,  convective  parameterization  to  prevent  a  supera- 
diabatic  layer  from  developing  at  the  nuxed  layer  top,  and  a  nudging  term  for  assimilating 
Doppler-radar-observed  horizontal  winds. 

In  order  to  run  the  ctunulus  parameterization  of  the  RAMS  model  (this  is  distinct 
from  the  convective  parameterization  for  the  mixed  layer  top),  3- dimensional  fields  are 
reconstructed  fixxn  the  prognostic  and  diagnostic  mixed  layer  fields  and  their  assumed 
vertical  structure,  and  are  written  into  the  regular  RAMS  model  3-dimensional  fields  in  the 
applicable  vertical  levels.  The  cumulus  parameterization  is  then  nm  directly  from  the  latter 
fields. 

The  mixed-layer  option  was  designed  to  be  totally  tmobtnisive  on  the  RAMS  model; 
it  requires  nothing  from  the  basic  model  which  is  not  already  provided  except  for  a  few 
additional  parameters  defined  in  the  input  namelist.  No  new  arrays  need  to  be  dimensioned 
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to  run  in  the  MLM  mode  as  sufficient  array  space  is  already  available.  Since  running  ILA.MS 
in  MLM  mode  turns  off  the  prognostic  equations  for  the  3-dimensional  model  variables,  the 
RAMS  tendency  curays  FU,  FV,  FW,  FP,  and  FTH  are  not  needed  elsewhere.  The  MLM 
code  has  thus  adopted  these  arrays  for  nearly  all  its  storage  requirements.  A  minimum  of 
six  vertical  levels  and  nine  gridpoints  in  each  horizontal  direction  mxist  be  dimensioned  in 
RAMS  to  provide  adequate  space  for  the  MLM. 

Many  of  the  regular  RAMS  namelist  variables  have  no  function  in  the  MLM  and  may 
be  ignored  when  running  in  that  mode.  Those  not  needed  in  the  MLM  are  indicated  in  the 
namelist  documentation. 

Some  comments  concerning  the  variables  not  needed  to  run  RAMS  in  MLM  mode  are 
the  following.  The  MLM  has  no  soil  model  or  radiation  parameterization,  and  uses  its 
own  time-  and  space-dependent  surface  duxes  of  heat  and  moisture  specified  by  the  user. 
Since  the  MLM  runs  many  times  faster  than  a  regular  RAMS  run  with  several  vertical 
levels,  it  was  decided  that  a  history  write  and  restart  would  probably  never  be  needed  so 
the  MLM  variables  have  not  been  set  up  to  be  written  to  a  history  file.  The  MLM  has 
no  small  timestep  caJculations  and  no  need  for  special  computation  to  minimize  latered 
boundary  reflection.  The  only  reason  DELTAZ,  DZRAT,  amd  DZMAX  are  needed  is  for 
providing  regular  3-dimensiontJ  RAMS  model  levels  to  be  filled  with  IvlLM  data  and  used 
for  diagnosing  cumulus  convection. 

3  Validation  simulations 

Following  the  construction  phase  of  the  two  model  versions  developed  for  this  project,  a 
series  of  tests  was  conducted  to  verify  that  each  version  performed  properly.  These  tests 
were  most  important  in  the  case  of  the  mixed-layer  model  since  it  was  started  from  scratch 
from  equations  and  other  descriptions  provided  in  the  literature,  and  grafted  into  the  RAMS 
code.  The  full  model  version,  on  the  other  hand,  was  constructed  as  a  customized  subset 
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of  FLA^S  which  is  a  tested  and  established  model.  Nevertheless,  verification  tests  were 
performed  on  the  full  model  as  well.  These  also  served  as  benchmarks  against  which  to 
compare  the  mixed-layer  model  solutions. 

Basic  simulations  with  the  mixed-layer  model  were  first  performed  to  test  every  individ¬ 
ual  physical  and  computational  process  represented  in  the  code.  The  processes  specifically 
performed  in  the  new  mixed-layer  model  include  the  following. 

1.  Reading  parameters  set  for  the  model  in  the  namelist  data 

2.  Constructing  from  the  RAMS  input  sounding  a  base  atmospheric  state  consistent  with 
the  imderlying  assumptions  on  the  vertical  structure  of  the  atmosphere  which  form 
the  basis  o(  the  mixed-layer  model 

3.  Initializing  all  prognostic  variables  used  in  the  mixed-layer  model  including  a  balanced 
boundary  layer  mean  F.lcmAn  profile  computed  iteratively. 

4.  Computing  advective  tendencies  for  all  prognostic  variables 

5.  Computing  the  contribution  to  the  mixed-layer  horizontal  pressure  gradient  force  due 
to  horizontal  variations  in  the  mixed-layer  potential  temperature,  the  inversion  top 
potential  temperatxire,  the  mixed-layer  top  height,  and  the  geostrop  hie  wind 

6.  Computing  the  Coriolis  force  tendency  on  the  horizontaJ  velocity  components 

7.  Computing  the  contributions  to  the  tendencies  of  mixed-layer  momentum,  moisture, 
and  potential  temperatme  due  to  surface  fluxes  of  each 

8.  Computing  the  contributions  to  the  tendencies  of  mixed-layer  momentum,  moisture, 
potential  temperature,  and  depth  due  to  entrainment  fliixes  at  the  top  of  the  mixed 
layer 

9.  Computing  horizontal  diffusive  tendencies  for  all  prognostic  variables 
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10.  Accumulating  all  prognostic  tendencies  and  advancing  prognostic  variables  forward  in 
time 

11.  Performing  convective  adjustment  on  prognostic  variables  any  time  the  top  of  the 
mixed  layer  reaches  a  superadiabatic  lapse  rate  (the  paramieterizations  in  the  model 
prevent  this  from  occurring  in  nearly  any  situation) 

12.  Diagnosing  vertical  velocity  and  other  quantities  from  the  prognostic  variables  and 
varioiu  model  specifications 

13.  Computing  appropriate  values  for  the  prognostic  variables  at  the  Isiteral  boimdaries 

14.  Filling  up  the  full  3-D  version  of  the  model  from  prognostic  variables  and  their  stssmned 
vertical  structure  in  the  mixed-layer  model  so  that  fields  may  be  plotted  tuid  possible 
convection  may  be  diagnosed 

Each  of  these  processes  was  carefully  checked,  using  the  VMS  debugger  and  sometimes 
a  hand  calculator  where  helpful.  Examples  of  the  simpler  tests  performed  are  the  following. 

1.  The  model  was  run  without  surface  fluxes  of  any  type  and  without  topography. 
Whether  the  model  was  initialized  with  wind  or  no  wind,  it  was  verified  that  the 
solution  remained  completely  balanced  and  constant  in  time. 

2.  The  model  was  nm  with  no  stirface  fluxes  except  for  momentiun.  The  input  sounding 
contidned  a  geostrophic  wind.  It  was  verified  that  the  balsmced  meeui  Ekman  wind  in 
the  mixed  layer  satisfied  the  equation  for  the  mean  Ekman  wind  for  various  surface 
drag  coefficient  values,  and  remained  balanced  and  constant  in  time  during  a  model 
integration. 

3.  The  model  was  initialized  with  no  wind.  A  surface  heat  flux  was  specified  and  the 
growth  of  the  mixed  layer  was  simulated.  It  was  verified  that  the  model-computed 
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entrainment  fluxes  at  the  top  of  the  mixed  layer  obeyed  the  referenced  equations  used 
to  design  the  model,  and  that  the  mixed  layer  growth  rate  and  temperature  increase 
were  consistent  with  the  surface  hesiting  rate. 

4.  The  evolution  of  the  mixed-layer  moisture  field  was  verified  for  consistency  with  the 
surface  and  top  entrainment  fluxes  of  moisture  in  a  no-wind  simulation. 

5.  The  individual  contributions  in  the  hydrostatically-computed  pressure  field  in  the 
mixed  layer  were  checked  at  individual  grid  points  against  the  original  equations  for 
a  model  state  containing  a  complicated  thermal  structtire. 

6.  Both  the  entrainment  and  convective  adjustments  made  to  the  prognostic  variables 
were  checked  to  verify  that  the  at^pxstments  were  being  made  correctly. 

Once  these  basic  checks  were  completed,  regular  simulations  of  fairly  standard  processes 
were  performed  to  see  if  the  expected  solution  would  develop.  We  report  on  two  such 
simulations  here.  The  first  was  a  sea  breeze  cirailation,  set  up  as  a  result  of  the  horizont2d 
temperature  contrast  between  a  constant  temperature  water  surface  tmd  a  rapidly  heating 
dry  land  surface  on  a  summer  day.  Because  such  a  circulation  will  tend  to  advance  toward 
the  land  surface,  it  is  customary  in  test  simulations  to  initialize  the  model  with  a  weak 
wind  blowing  Cram  land  to  sea.  This  was  done  here  using  a  2  m/s  wind  blowing  from  the 
east,  the  land  occupying  the  eastern  half  of  the  domain  and  the  sea  the  western  half.  The 
computational  domain  used  in  both  models  spanned  200  km  in  the  east- west  direction  and 
100  km  in  the  north-south  direction  with  5  km  resolution  in  both  directions. 

The  initial  potential  temperature  field  used  in  both  models  is  shown  in  a  vertical/east- 
west  cross  section  through  the  center  of  the  domain  in  Figure  3.1.  It  is  weakly  stable  up  to 
a  height  of  3  km  and  strongly  stable  above  that.  After  a  simulation  period  of  4  hours,  the 
land  surface  has  warmed  considerably  while  the  sea  surface  has  not.  This  is  illustrated  in 
the  potential  temperature  fields  for  both  the  full  and  mixed-layer  versions  of  the  model  in 
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figures  3.2a  and  3.3a,  respectively.  The  surface  heat  fluxes  in  the  two  models  were  not  made 
equal,  so  the  full  model  has  undergone  more  land  warming.  The  highly  baroclinic  region 
near  x  =  0  dnves  a  solenoidal  circiilation  illiistrated  in  the  u-component  and  w-component 
velocity  fields  in  figures  3.2b  and  3.2c  for  the  full  model  and  in  figures  3.3b  and  3.3c  for 
the  mixed-layer  model.  The  circulation  is  considerably  stronger  in  the  full  model  because 
of  the  stronger  land  heating,  but  both  models  exhibit  a  similetf  structure  in  the  horizontad 
wind  field.  Note  that  the  mixed-layer  model  does  not  have  an  explicit  return  flow  above 
the  mixed-layer,  so  only  the  single  layer  wind  is  shown.  The  v-component  (the  north-south 
direction,  into  the  page)  wind  fields  for  the  two  models  are  shown  in  figures  3.2d  and  3.3d. 
The  v-component  has  developed  under  the  influence  of  the  Coriolis  force  acting  on  the 
solenoidal  circtilation.  The  main  difference  between  the  two  cases  for  this  field  is  that  the 
mixed-layer  solution  lacks  the  positive  v-component  region  at  the  two-kilometer  level  just 
west  of  the  domain  center.  This  is  another  consequence  of  the  absence  of  an  explicit  return 
flow  above  the  mixed  layer. 

The  second  regular  simulation  reported  here  contained  a  solenoidal  circulation  driven  by 
strong  surface  heating  of  sloping  terrain.  This  simulation  used  a  150  km  by  150  km  horizon¬ 
tal  domain  in  which  the  southwest  fourth  of  the  domain  contained  a  mesa-like  topographic 
formation  having  a  horizontal  top  and  straight  sloping  surfaces  on  its  north  and  east  sides. 
A  2.8  m/s  southwesterly  flow  was  specified  as  the  initial  soimding  wind,  and  the  vertical 
temperature  profile  was  initialized  as  in  the  sea  breeze  experiments  above.  Figures  3.4a 
and  3.4b  show  a  hoiizontal  cross  section  of  the  mixed-layer  horizontsJ  winds  two  hours  into 
the  simulation.  The  u-component,  which  initially  was  2.0  m/s  throughout  the  domzun  has 
decreased  everywhere,  attaining  a  value  of  -2.0  m/s  along  the  eastern  side  of  the  mesa.  This 
indicates  a  well-developed  upslope  flow.  The  v-component  shows  a  similar  pattern  on  the 
north  slope  of  the  mesa.  However  it  only  achieves  a  minimum  value  near  0  m/s  and  actuedly 
increases  above  its  initial  2.0  m/s  value  in  the  northern  and  eastern  parts  of  the  domain. 
The  asymmetry  between  u  and  v  is  due  to  the  Coriolis  effect.  The  southwestern  location  of 
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the  higher  terrain  causes  the  mean  upslope  circulation  to  flow  toward  the  southwest.  The 
Coriolis  force  acting  on  this  ageostrophic  perturbation  from  the  initial  flow  is  toward  the 
northwest.  Hence,  the  Coriolis  force  acts  to  increase  the  v-component  and  decrease  the  u- 
component.  The  mesa  rim,  being  the  location  where  the  horizontal  Lapladan  of  the  terrain 
height  is  negative,  is  a  region  of  surface  convergence  and  upward  vertical  motion.  This  can 
be  seen  in  the  w-component  of  motion  shown  in  flgure  3.4c.  At  the  northeast  comer  of  the 
mesa  top,  where  the  northern  and  eastern  rims  meet,  the  upward  motion  is  enhanced.  All 
these  results  are  physically  expected. 

4  Testing  of  wind  data  assimilation 

A  primary  goal  of  the  software  development  phase  of  this  project  was  to  design  and  im¬ 
plement  a  technique  for  assimilating  wind  data  obtained  firom  single  Doppler  radar  into 
the  numerical  prediction  models.  After  consulting  various  references  on  the  subject,  we 
decided  that  a  nudging  method  would  be  appropriate  for  both  versions  o(  the  model.  The 
nudging  algorithm  is  implemented  by  adding  a  forcing  term  to  both  horizontal  equations 
of  motion.  This  term  appears  for  each  component  in  the  form  d/dt(VMOD)  =  (VOBS- 
VMOD)/SCLNUD  where  VOBS  is  the  observed  value  of  the  component  obtained  from  the 
Doppler  radar,  VMOD  is  the  current  model  vzJue  of  the  component,  and  SCLNTJD  is  a 
time  scale  determining  the  rate  at  which  VMOD  is  to  be  forced  toward  VOBS.  Typically, 
SCLNUD  would  be  on  the  order  of  1/2  hour,  a  time  short  enough  for  relaxation  to  take 
place  in  a  fraction  of  the  total  simulation  time  and  long  enough  so  that  the  model  is  not 
completely  forced  to  the  observed  wind  but  is  allowed  to  retain  some  of  its  own  history. 
We  developed  nudging  algorithms  for  assimilating  both  cartesian  wind  components  derived 
from  the  single  Doppler  analysis  technique  described  by  Easterbrook  (1975)  and  radial 
winds  obtained  directly  from  the  Doppler  measurements.  The  tests  of  the  first  algorithm 
are  described  in  this  chapter,  and  the  tests  of  the  second  in  Chapter  5. 
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The  primary  intended  application  of  the  project  under  which  these  models  were  devel¬ 
oped  is  forecasting  the  first  occurrence  of  significant  convection  within  a  specified  domain. 
The  forecast  is  to  be  aided  by  horizontal  radial  wind  measurements  made  by  a  single  Doppler 
radar  which  is  scanning  the  planetary  boundary  layer  in  the  domadn.  The  most  common 
boundary  layer  wind  feature  which  would  accompany  the  early  development  of  strong  con¬ 
vection  is  horizontal  convergence  with  vertical  divergence.  We  therefore  concentrated  on 
testing  the  assimilation  of  wind  data  into  the  model  for  cases  containing  a  region  of  hori¬ 
zontal  convergence  in  the  wind  field.  The  most  straightforward  way  of  defining  such  a  wind 
field  was  to  specify  a  convergence  line  parallel  to  one  of  the  coordinate  axes  (we  chose  a 
north-south  orientation)  through  the  center  of  the  domain.  Only  the  u-component  of  mo¬ 
tion  (that  perpendicular  to  the  convergence  line)  was  assumed  to  be  nonzero,  and  this  was 
defined  to  be  spatially  constant  on  either  side  of  the  line  while  varying  linearly  across  the 
line. 

One  such  test  performed  with  both  versions  of  the  model  developed  for  this  project  is 
presented  here.  The  models  were  initialized  with  a  sotmding  representing  a  relatively  stable 
atmosphere  with  no  wind.  The  model  domains  were  each  200  km  in  the  east-west  direction 
and  100  km  in  the  north-south  direction  with  a  5  km  grid  spacing  in  both  directions.  The 
convergence  line  in  the  assimilated  data  was  set  to  a  width  of  10  km  or  two  grid  points, 
with  3  m/s  westerly  winds  to  its  west  and  3  m/s  easterly  winds  to  its  east.  The  nudging 
relaxation  time  was  set  to  1800  s  in  both  models. 

Figures  4.1  and  4.2  show  the  model  solutions  for  the  full  and  mixed-layer  versions  of 
the  model  one  hour  into  the  simulation.  Figures  4.1a  and  4.2a  present  a  verticaJ/east-west 
cross  section  through  the  middle  of  the  domain  of  the  u-component  of  motion.  Both  cases 
show  a  definite  region  of  surface  convergence  centered  at  the  x  =  0  location,  but  the  width 
of  the  region  is  far  greater  than  the  10  km  width  of  the  assimilated  data.  Note  that  the  full 
model  depicted  in  figure  4.1  contains  a  compensating  horizontal  divergence  aloft  while  the 


13 


mixed-layer  model,  which  has  no  explicit  return  flow,  shows  only  the  convergent  mixed-layer 
winds.  Figures  4.1b  and  4.2b  show  the  v-component  which  has  developed  as  a  result  of  the 
Coriolis  force  acting  on  the  u-component,  which  is  ageoetrophic  because  it  was  not  present 
in  the  initial  soiinding.  The  vertical  motion  fields  which  are  consistent  with  the  horizontal 
convergence  fields  in  the  two  models  are  shown  in  figures  4.1c  and  4.2c. 

The  fact  that  the  u-component  has  not  developed  a  narrow  convergence  feature  like  the 
assimilated  data,  even  though  the  simulation  has  progressed  for  two  nudging  relaxation  time 
scale  periods,  implies  that  some  other  forcing  term  in  the  u-component  eqiiation  of  motion 
is  comparable  in  magnitude  to  the  nudging  term.  We  found  this  to  be  the  pressure  gradient 
force  term.  The  perturbation  pressiire  field  is  shown  in  figure  4. Id  for  the  full  model.  The 
pressure  has  risen  a  few  millibars  in  the  domain,  most  strongly  at  the  surface.  The  presstire 
rise  is  the  direct  result  of  adiabatic  cooling  of  the  stable  air  as  it  is  forced  upward  in  response 
to  the  nud^ng  the  horizontal  wind  toward  convergence.  The  cooling  is  most  prominent 
in  the  center  of  the  domain  where  upward  motion  is  strongest,  as  seen  from  a  plot  of  the 
potential  temperature  fields  in  figures  4.1e  and  4. 2d.  As  a  result,  the  surface  pressure  is 
highest  at  that  location  as  well.  The  pressure  gradient  force  is  thus  directed  away  from  the 
domain  center  on  both  sides  which  opposes  the  nudging  force.  Hence,  the  sharp  definition 
of  the  convergence  line  in  the  assimilated  data  is  prevented  from  developing. 

Four  hours  into  the  simulation,  the  difference  between  the  model  winds  and  the  usimi- 
lated  winds  is  even  larger.  Figures  4.3  and  4.4  show  the  full  model  and  mixed-layer  model 
results,  respectively,  for  this  later  time.  The  u-component  fields  of  figures  4.3a  and  4.4a 
show  an  even  broader  zone  of  convergence,  with  the  u-component  varying  by  less  than  5 
m/s  over  a  200  km  distance.  This  is  quite  different  from  the  6  m/s  change  over  a  10  km 
distance  present  in  the  assimilated  u-component  field.  The  full  model  perturbation  pressure 
field  plotted  in  figure  4.3c  shows  the  horizontal  pressure  gradient  force  at  the  surface  to  be 
stronger  than  before.  The  vertical  velocity  fields  for  both  models  plotted  in  figures  4.3b  and 
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4.4b  show  the  region  of  upward  motion  to  have  expanded  significantly  in  the  x-direction. 
The  strongest  upward  motion  still  remains  in  the  center,  however.  A  corresponding  major 
lateral  expansion  of  the  cold  dome  in  both  models  is  likewise  apparent  from  the  potentied 
temperature  fields  in  figures  4.3d  and  4.4c. 

The  problem  worsens  when  the  simtilation  is  carried  out  further  with  the  nudging  turned 
off.  Because  the  intended  application  of  these  models  is  to  mcike  a  forecast,  normal  use  of 
them  should  involve  sissimilating  wind  data  for  a  period  of  time  and  then  continuing  the 
simulation  in  pure  forecast  mode  without  continuing  the  data  assimilation.  This  was  done 
for  the  present  test  cases  by  tximing  off  the  nudging  immediately  afler  the  four  hour  time, 
amd  continuing  the  simulations  for  another  hour.  Figures  4.5a  and  4.6a  show  the  full  model 
and  mixed-layer  model  u-component  fields  after  the  additional  hour.  Instead  of  surface 
convergence,  a  significant  region  of  divergence  has  formed  in  the  middle  of  the  domain. 
The  corresponding  vertical  velocity  fields  shown  in  figures  4.5b  and  4.6b  show  subsidence  to 
pervade  the  middle  third  of  the  domcun.  The  reason  for  this  is  that  the  cold  dome  built  up 
under  the  influence  of  the  nudging  force  is  now  in  the  process  of  collapsing.  A  comparison 
of  the  full  model  theta  field  shown  for  this  time  in  figure  4.5c  with  the  earlier  one  in  figure 
4.3d  shows  the  central  cold  peak  to  have  disappeared.  That  the  surface  winds  are  divergent 
rather  than  convergent  at  this  time  shows  that  the  nudging  of  the  model  to  the  observed 
wind  fields  has  proved  a  detriment  rather  than  an  improvement  to  the  solution. 

In  summary,  the  models  have  failed  to  respond  to  the  nudging  force  in  such  a  way  that  a 
new  bfdance  of  the  model  fields  is  achieved  and  the  model  velocity  field  agrees  closely  with 
the  assimilated  data.  This  failure  is  a  consequence  of  the  assimilated  data  containing  wind 
data  inconsistent  with  the  models’  temperature  fields  or  mass  fields.  In  a  stable  atmosphere, 
as  used  for  these  tests,  the  temperature  and  velocity  fields  are  intimately  coupled  so  that 
a  change  in  one  induces  a  change  in  the  other.  In  turn,  the  induced  change  in  the  second 
induces  a  coimter  change  in  the  first,  opposing  the  original  change  to  the  first. 
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If,  on  the  other  hand,  a  nearly  neutral  atmosphere  is  used  to  initialize  the  models,  the 
temperature  field  is  loosely  coupled  to  the  velocity  field  so  that  any  changes  induced  in  the 
velocity  field  do  not  produce  changes  in  the  mass  field  which  fight  the  velocity  chamges.  We 
performed  simulations  with  both  models  which  demonstrate  this  point.  These  simulations 
are  identical  to  the  stable  ones  described  above,  except  that  they  were  initialized  with  neairly 
neutral  lapse  rates  up  to  a  height  o{  3  km.  Figures  4.7  and  4.8  show  the  u-component  fields 
after  4  hours  of  simulation  (and  nudging)  time.  In  this  case,  both  models  have  been  able 
to  develop  and  maintain  fairly  strong  convergence  lines,  although  they  are  stiU  not  as 
sharp  as  in  the  assimilated  data.  As  before,  the  nudging  is  turned  off  at  the  four  hour 
time  and  the  simtdation  continued  for  one  additiortal  hour.  The  full  model  u-component 
field  (figure  4.9)  still  maintains  surface  convergence,  although  it  has  weadtened.  The  mixed 
layer  model  u-component  field  (figure  4.10)  however  is  exhibiting  some  surface  divergence. 
This  is  a  consequence  of  the  model's  lack  of  an  explicit  return  flow  above  the  mixed-layer 
convergence  field.  The  mixed-layer  model  makes  more  approximations  than  the  full  model 
and  hence  introduces  more  adjustable  parameters.  Appropriate  tuning  of  these  parameters 
is  required  for  the  model  response  to  more  closely  follow  the  behavior  of  the  full  model.  In 
this  case,  using  a  larger  value  oi  the  horizontal  diffusion  coefficient  for  the  boundary  layer 
height  would  have  prevented  it  from  growing  so  rapidly  in  the  presence  of  the  sustained 
vertical  motion  in  the  center  of  the  domain.  Because  the  model  implicitly  requires  the 
mixed  layer  to  be  capped  by  a  stable  layer,  the  higher  diffusion  would  allow  the  upweird 
motion  to  continue  without  a  continued  elevating  of  the  mixed  layer  top  and  a  consequent 
building  of  a  cold  dome. 

The  main  finding  of  these  generic  data  assimilation  tests  is  that  as  long  as  the  tem- 
perattire  and  velocity  fields  are  coupled  due  to  atmospheric  stability,  it  is  not  sufficient  to 
assimilate  only  velocity  data  into  a  numerical  prediction  model;  appropriate  temperature 
data  are  also  required.  Although  a  negative  result,  it  is  an  important  one  in  view  of  the 
fact  that  there  are  widespread  plans  to  assimilate  the  velocity  fields  derived  from  the  up- 
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coming  national  network  of  NEXRAD  radars  and  wind  profilers  into  operational  forecast 
models  without  temperature  information  of  the  same  detail  and  coverage.  This  fact  amd 
the  negative  results  we  obtained  in  the  nudging  tests  for  this  study  prompted  us  to  submit 
a  paper  to  the  Radar  Conference  which  is  to  be  held  in  March  1989.  In  the  paper,  we 
describe  a  physical  scenario  in  which  a  single  Doppler  radar  might  detect  a  convergence  line 
in  the  wind  field  but  not  be  able  to  ascertain  the  thermal  structure  of  the  line.  A  simulation 
to  demonstrate  the  points  made  above  was  performed  using  the  nonhydrostatic  version  of 
RAMS  which  executes  more  slowly  than  the  models  described  here  but  provides  full  detail 
and  directly  interpretable  results.  The  paper  is  included  as  Appendix  B  of  this  report.  We 
conclude  that  to  effectively  assimilate  into  a  numerical  simulation  observed  wind  data  which 
implies  vertical  motion,  it  is  necessary  to  simultaneously  assimilate  temperature  data  which 
is  consistent  with  that  vertical  motion  field. 

5  Testing  of  radial  nudging  algorithm 

We  developed  and  tested  for  this  project  a  new  algorithm  for  assimilating  single  Doppler 
wind  data  based  on  a  concept  proposed  by  J.  Goerss,  formerly  of  the  Cooperative  Institute 
for  Mesoscale  Metecnological  Studies,  emd  R.  Doviak  emd  R.  Rabin  of  the  National  Severe 
Storms  Laboratory  in  Norman,  Oklahoma.  Their  concept  was  to  directly  input  to  the 
model  the  observed  single  Doppler  radial  wind  component,  rather  than  first  processing  the 
observed  radial  wind  through  another  technique  (e.g.,  Easterbrook,  1975)  to  extract  the 
two  horisontal  cartesian  wind  components. 

Obviously,  any  technique  for  computing  the  totsJ  horizontal  wind  field  from  a  radial 
component  only  is  adding  information  on  the  flow.  Some  of  this  information  is  accurate 
because  it  is  derived  from  constraints  imposed  on  the  flow  by  the  basic  laws  of  fluid  me¬ 
chanics  such  as  mass  continuity.  However,  these  laws  are  not  sufficient  to  close  the  set  of 
equations  used  in  obtaining  the  cartesian  components,  and  additionsd  arbitrary  assump- 
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tions  are  required.  As  a  consequence,  significant  imcertainty  exists  in  the  final  product. 
By  eliminating  this  step  in  the  transfer  of  Doppler  information  to  the  model  and  directly 
inputing  the  radial  component,  the  introduction  of  extraneous  information  to  the  winds  is 
also  eliminated. 

Because  the  model  employs  the  equiUions  of  fluid  dynamics,  it  is  hypothesized  that  it 
will  place  appropriate  physical  constriunts  on  the  assimilated  radial  winds  and  effectively 
extract  from  them  both  cartesian  components.  This  technique  has  a  chance  of  being  superior 
to  the  standard  method  above  in  that  the  missing  information  needed  in  obtaining  the  full 
horizontal  wind  field  from  the  radial  component  is  in  effect  supplied  by  the  model  wind  itself, 
and  is  thus  not  arbitrary  but  is  instead  related  to  the  physical  situation  being  observed. 

We  designed  a  numericsd  algorithm  which  nudges  the  model  winds  directly  to  the  input 
radial  component.  The  method  consists  of  the  following  steps. 

1.  At  each  grid  point  in  the  model  domain,  the  current  model  cartesian  components 
are  decomposed  into  a  radial  and  azimuthal  component  based  on  a  polar  coordinate 
system  centered  at  the  location  of  the  Doppler  radar. 

2.  The  difference  between  the  observed  radial  wind  and  the  model  radial  wind  is  com¬ 
puted,  and  the  model  rstdial  wind  is  nudged  toward  the  observed  radicd  wind  with 
a  strength  proportional  to  that  difference.  The  model  azimuthal  component  is  not 

altered. 

3.  After  the  forcing  of  the  model  radial  component,  the  radial  and  azimuthal  components 
are  transformed  back  into  the  cartesian  components. 

We  tested  this  radial  nudging  algorithm  for  several  wind  fields  in  which  both  cautesiain 
components  were  specified  according  to  some  idealized  flow  pattern.  The  r2Mlial  wind  com¬ 
ponent  which  a  Doppler  radar  would  measure  from  each  of  these  fields  was  computed  and 
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used  as  the  observed  wind  for  assimilation  into  the  model.  The  model  winds  were  then 
examined  throughout  their  evolution  to  determine  how  closely  the  radial  nudging  algorithm 
forced  them  to  the  initial  cartesian  wind  components. 

The  first  tests  we  conducted  of  this  algorithm  involved  wind  fields  containing  horizontal 
convergence  somewhere  in  the  domain.  These  produced  quite  unacceptable  results  auid  led 
us  to  suspect  that  the  method  was  deficient  and  would  not  work.  We  later  discovered 
that  even  straightforward  cartesian  nudpng  of  both  wind  components  will  fail  when  it  tries 
to  produce  significant  vertical  motion  in  stable  air  and  thermodynamic  data  which  would 
conform  to  this  vertical  motion  is  not  available  for  simultaneous  assimilation  (see  Section 

4  and  Appendix  B  of  this  report).  We  thus  re-tested  the  radial  nudging  algorithm  on  a 
horizontal  wind  field  containing  translation,  rotation,  and  deformation,  but  no  divergence. 
Moreover,  we  input  a  stable  atmosphere  for  the  test  to  effectively  constrain  the  horizontal 
wind  field  to  remain  nondivergent.  This  constraint  increases  the  probability  of  success  of  the 
algorithm  for  this  case  because  the  correct  solution  is  itself  nondivergent.  The  radial  wind 
component  is  divergent,  even  though  the  cartesism  field  from  which  it  is  obtained  is  not,  and 
nudging  to  the  radial  component  tries  to  force  the  model  winds  to  become  divergent.  The 
large  atmospheric  stability  specified  for  this  case  prevents  the  vertical  motion  and  hence 
the  horizontal  divergence  from  ever  becoming  large  for  a  sustained  period  of  time. 

The  radial  nudpng  test  simulation  was  run  using  the  mixed-layer  version  of  RAMS  on 
a  physical  domain  spanning  200  km  in  the  x-direction  and  100  km  in  the  y-direction  with 

5  km  resolution  in  both  directions.  The  assxuned  location  of  the  Doppler  radar  is  very  near 
the  center  of  the  domain.  The  prescribed  nondivergent  wind  field  whose  radial  component 
will  be  assimilated  into  the  model  consists  of  a  uniform  flow  from  the  south  of  3  m/s  (v  = 
3  m/s)  over  the  eastern  two-thirds  of  the  domain  and  a  imiform  flow  from  the  north  of  3 
m/s  (v  =  -3  m/s)  over  the  western  one-third  of  the  domain.  The  u-component  of  motion 
is  zero  throughout  the  prescribed  wind  field.  The  two  constant  v  regions  are  separated  by 
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a  10-km-wide  shear  line  over  which  the  wind  varies  linearly.  The  entire  shear  line  contains 
vorticity  and  deformation,  but  no  mean  translation  and  no  divergence.  The  regions  on 
either  side  of  the  shear  line  contain  pure  translation  only.  Nudging  of  the  model  to  the 
radial  wind  component  of  this  field  was  performed  over  a  four-hour  simulation  period  using 
a  nudging  relaxation  time  scale  of  1800  s.  Nudging  was  then  ceased  and  the  simulation 
carried  out  for  an  additional  two  hours  to  examine  any  flow  evolution  which  might  occur  in 
the  absence  of  the  forcing. 

Figures  S.la-c  show  the  u,  v,  and  w  components  of  motion  1  hour  after  the  simulation 
(and  nudging)  startup  time.  Any  nonxero  value  for  the  u-component  is  of  course  a  deviation 
from  the  actual  cartesian  wind  field  from  which  the  assimilated  radial  component  is  com¬ 
puted.  It  arises,  however,  as  a  consequence  of  the  radial  nudging  algorithm  which  has  the 
requirement  that  nud^ng  wUl  not  directly  force  the  azimuthal  component.  As  long  as  the 
v-component  is  to  be  directly  forced  by  the  nudging  at  a  location  not  due  north  or  south  of 
the  radar,  the  u-component  must  be  fenced  simultaneously  to  satisfy  that  reqtiirement.  Note 
that  the  u-component  northwest  of  the  radar  location  (which  is  near  x=0,y=0)  is  negative 
implying  flow  toward  the  west.  The  v-component  is  positive  in  the  same  location  so  that 
the  resultant  flow  there  is  radially  outward  from  the  radar.  This  is  in  direct  accordance 
with  the  outward  radial  motion  that  the  radar  would  see  in  the  total  cartesian  wind  field 
specified  for  this  experiment. 

A  wmilay  discussiou  applies  to  the  other  regions  of  the  flow  in  figures  5.1a  and  5.1b.  This 
tim*  in  the  simulation  is  rather  early  for  horizontal  pressure  gradient  forces  to  have  exerted 
a  significant  influence  on  the  horizontal  velocity  field,  so  the  latter  has  thus  far  responded 
primarily  to  the  nudging  forces  only.  The  horizontal  flow  does  however  contain  substantial 
rates  of  divergence  and  convergence  as  deducible  from  figure  5.1c.  Vertical  velocities  have 
exceeded  2  cm/s  in  magnitude  at  a  height  of  400  m  by  this  time.  The  location  of  strongest 
convergence  can  be  seen  to  be  due  south  of  the  radar.  This  is  producing  surface  high  pressure 
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due  to  the  stability  of  the  air  and  its  ascent  above  the  convergence  which  will  counteract 
the  positive  u-component  to  its  west  (figure  5.1a).  Similjuly,  a  surface  low  pressure  is 
forming  due  north  of  the  radar  because  of  the  subsidence  wanning  above  that  location,  the 
the  negative  u-component  seen  to  its  west  in  figure  5.1a  will  be  coimteracted  as  a  result. 
Farther  west  of  these  locations,  just  west  of  the  shear  line  at  x  =  -40  km,  the  north-south 
pattern  of  upward  and  downward  motion  is  reversed,  as  is  the  pattern  of  the  positive  and 
negative  u-component  motion,  so  there  the  developing  horizontal  pressure  gradient  forces 
will  again  tend  to  force  the  u-component  in  the  direction  of  zero.  The  v-component  in  figure 
5.1b  is  showing  a  clear  development  of  the  shear  zone  northwest  and  southwest  of  the  radar 
where  it  is  directly  detected  in  the  radial  wind  component,  but  not  west  of  the  radar  where 
the  radial  component  contains  no  information  on  the  zone  at  all.  However,  the  north-south 
component  of  the  pressure  gradient  force  which  is  developing  on  each  side  of  the  shear  zone 
(in  the  opposite  sense  on  each  side)  wiU  strengthen  the  shear  west  of  the  radar  in  its  forcing 
of  the  v-component.  This  demonstrates  how  the  constraint  imposed  on  the  flow  by  the  large 
static  stability  which  was  specified  for  this  case  is  helping  to  force  the  model  flow  to  the 
true  solution. 

Figures  5.2a-c  show  the  u,  v,  and  w  components  of  motion  in  the  model  at  fotir  hours 
after  the  startup  time  of  the  simulation.  Following  the  initial  transient  induced  by  the 
midgjtifl  early  in  the  simulation,  the  verticsJ  velocity  field  has  damped  to  weaker  values, 
not  exceeding  0.5  cm/s  in  most  of  the  domsdn.  All  three  components  have  become  fairly 
steady,  not  having  changed  much  from  the  3-bour  simulation  time  (not  shown).  Thus,  the 
model  has  sg)prosiched  what  will  be  its  final  solution.  The  u-component  in  this  solution 
has  a  systematic  tendency  toward  eastward  motion  in  the  north  half  the  domain  and 
westward  motion  in  the  south  half  giving  a  negative  mean  vorticity  (for  that  component) 
which  is  opposite  to  the  sense  of  the  specified  shear  zone  vorticity.  The  v-component,  on 
the  other  hand,  exhibits  a  mean  positive  vorticity  over  the  domain,  the  same  sense  as  that 
specified  in  the  shear  zone.  Compared  with  figure  5.1b,  figure  5.2b  shows  the  shear  zone 
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to  be  continuous  from  north  to  south  without  a  breedc  west  of  the  radar.  The  filling  in 
of  the  middle  region  was  caused  by  the  pressure  gradient  forces  as  disoissed  above.  The 
v-component  overall  resembles  the  true  solution  fairly  well,  although  the  east  and  west 
extremes  of  the  domain  both  contain  erroneously  wetik  wind  magnitudes. 

Figures  5.3a-c  show  the  u,  v,  and  w  components  of  motion  at  the  six-hour  time  in  the 
simulation,  which  is  following  two  hours  without  nudging.  The  u-component  continues  to 
exhibit  a  negative  contribution  to  the  mean  vorticity  in  the  domain  while  the  v-component 
continues  to  do  the  opposite.  The  v-component  has  retained  a  well-defined  shear  zone 
without  radial  wind  data  assimilation,  although  the  zone  is  more  diffuse  across  its  width 
and  the  shear  line  has  begun  to  take  on  a  south- southwest  to  north-northeast  orientation 
in  response  to  advection  by  the  u-component.  The  w-component  has  undergone  the  largest 
relative  change  in  the  two-hour  period,  but  magnitudes  remain  small,  averaging  around  2 
or  3  mm/s  in  absolute  value.  The  horizontal  components  have  changed  only  slightly  during 
the  two  hours  since  the  fields  shown  in  figtires  5.2a-b  indicating  that  by  the  fom-hour  time, 
the  nudging  force  was  no  longer  strong  (i.e.,  that  the  model  radial  wind  field  had  become 
close  to  the  observed  radial  wind).  Thus,  the  model  attained  a  nearly  steady  stsite  solution 
different  from  the  specified  true  sc^ution. 

We  give  the  radial  nudging  simulation  a  fair  rating  for  this  case,  but  we  reiterate  that 
the  physical  scenario  was  deliberately  chosen  to  increase  the  likelihood  of  success.  We  have 
seen  that  the  technique  does  not  perform  as  well  in  the  general  case,  primarily  because 
of  the  problems  encoxmtered  when  nudging  attempts  to  force  vertical  motion  in  a  stable 
atmosphere.  That  is  a  separate  problem  which  is  encountered  in  the  assimilation  of  both 
cvtesian  wind  components  as  well.  We  are  still  of  the  opinion  that  some  of  the  reasoning 
behind  the  radial  nudging  technique  is  sound.  In  particular,  the  introduction  of  erroneous 
information  to  the  wind  field  which  occurs  in  standard  single  Doppler  analysis  techniques 
is  an  ongoing  problem  to  which  a  solution  should  be  found.  We  have  shown  that  direct 


22 


assimilation  of  the  observed  radial  component,  which  bypasses  the  Doppler  amalysis  step, 
is  not  a  viable  £dtemative.  During  this  experimentation,  we  have  identified  a  n^w  method 
which  has  the  potential  of  solving  both  problems.  Both  techniques  discmsed  above  (prior 
recovery  of  cartesian  components  before  assimilation  and  direct  radial  wind  assimilation) 
inv<Jve  only  a  one-way  communication  of  information;  Doppler  data  is  input  to  the  model 
but  model  information  is  not  fed  back  to  the  Doppler  data  external  to  the  model.  A  possible 
new  method  would  be  to  use  the  current  model  winds  as  input  to  a  st£mdard  single  Doppler 
analysis  technique  in  place  of  making  any  arbitrary  assumptions  on  the  wind  field.  We  plan 
to  test  this  idea  in  future  research.  A  solution  stiU  needs  to  be  foimd  to  the  more  serious 
problem  that  assimilation  of  wind  data  without  thermodynamic  information  is  generally 
unsuccessful. 

6  Simulations  of  actual  atmospheric  cases 

Although  it  was  entirely  expected  that  the  results  shown  by  the  “generic”  nudging  simu¬ 
lations  described  in  chapter  4  of  this  report  would  be  applicable  to  a  simulations  run  with 
actual  observations,  the  two  versions  of  the  model  were  nm  for  three  case  studies  from  the 
CINDE  field  program.  The  following  section  will  describe  the  radar  analysis  technique  amd 
the  results  of  those  simulations. 

A.  Methodology  for  Analysis  of  Doppler  Radar  Data  for  Model  Input 

For  the  real  data  case  studies,  the  single  Doppler  radial  velocity  data  from  NCAR's 
5-cm  CP-3  radar  during  1987  CINDE  experiment  was  used.  Three  separate  days  were 
considered,  17  July,  27  June,  and  31  July.  For  the  first  case,  four  elevations  (0.4,  0.9, 
1.7  and  2.6  deg)  were  used,  while  only  the  0.4  deg  elevation  was  used  in  the  latter  two 
cases.  Editing  consisted  of  thresholding  out  low  signal- to- noise  returns  and  those  velocities 
with  excessive  spectral  width,  unfolding  velocities,  removing  ground  clutter  and  second  trip 
echoes,  and  applying  a  liberal  spatial  filter  to  the  data.  All  this  was  done  in  ‘radar  space 
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using  NCAR’s  RDSS  software. 


The  edited  radial  velocity  data  were  then  interpolated  to  a  cartesian  grid  using  NCAR’s 
REORDER  interpolation  software.  Since  CP-3’s  western  view  was  severely  limited  by  the 
FVont  Rtuige  at  about  30  km  to  the  west,  the  100  by  100  km  grid  was  configured  off-center 
from  the  radar  in  the  E-W  direction,  from  30  km  to  the  west,  to  70  km  to  the  east  (about  the 
limiting  range  of  boimdary-layer  detectability).  The  grid  extended  50  km  from  the  radar  in 
both  the  north  and  south  directions.  The  horizontal  grid  spacing  of  the  Doppler  data  grid 
was  2.5  km  to  account  for  the  ntimerical  model  grid  stagger.  The  winds  were  interpolated 
to  vertical  grid  levels  at  0.5,  1.0,  1.5  and  2.0  km  above  the  radar  elevation.  The  radar 
sampling  capabilities  and  interpolation  method  yielded  quite  detailed  and  significant  radial 
velocity  gradients  in  the  horizontal,  with  vertical  shear  much  less  resolved.  (Since  only  one 
elevation  was  used  in  the  latter  two  cases,  the  radial  wind  fields  were  identical  at  all  levels.) 
These  gridded  radial  velocities  were  filtered  further,  and  echo-free  patches  in  the  domain 
were  filled  by  horizontal  interpolation  and  extrapolation  to  the  domain  boxmds. 

FVom  the  radial  velocity  field,  a  divergence  field  could  be  diagnosed  by  the  method  of 
Koscielny  et  al.  (1982).  However,  the  models  were  not  designed  to  be  nudged  to  accordance 
to  the  observed  radial  wind  and/or  divergence  fields,  but  instead  required  the  Cartesian 
wind  component  fields.  While  we  developed  and  initially  tested  a  method  for  deriving  the 
ttmgential  wind  field  as  a  residual  of  the  radial  wind  and  divergence  fields  (which  would 
have  provided  the  u  and  v  wind  fields),  it  could  not  be  applied  to  these  cases  because  of  the 
off-centered  position  of  the  rtular  in  the  100  by  100  km  computational  domain.  Thus,  for 
simplification,  the  tangential  wind  field  was  derived  from  the  NCAR  PAM  mesonet,  which 
had  a  station  spacing  of  about  10  km.  While  the  radial  wind  field  varied  with  height  (for 
the  first  case  only),  the  tangential  wind  field  was  constant  with  height.  We  did  not  explore 
methods  by  which  the  tangential  wind  could  have  varied  with  height  according  to  Eckman 
turning.  FVom  these  radar-derived  radial  winds  and  PAM-derived  tangential  winds,  the  u 
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and  V  components  were  obtained  for  input  to  the  model. 


Figures  6.1a-f  illustrate  the  processing  methodology.  The  PAM  observations  for  the 
first  case  at  1500  local  time  (LT)  are  shown  in  Figure  6.1a,  along  with  the  100  by  100  km 
computational  domain  and  the  off-centered  CP-3  radar  position.  The  dominauit  feature  is 
the  north-south  convergence  line  about  30  km  east  of  the  radair,  denoted  by  dashed  line 
A-B.  Figure  6.1b  shows  the  PAM  wind-field  as  objectively  analyzed  onto  the  full  grid-point 
domain,  and  the  corresponding  PAM-derived  divergence  field  shows  the  strong  convergence 
line.  Figure  6.1c  shows  the  radar-derived  radial  wind  field  (for  the  lowest  vertical  level), 
and  Figure  6. Id  shows  the  contoured  tangential  component  of  the  PAM  wind  field.  These 
radial  and  tangential  components  are  combined  to  yield  the  final  wind- vector  field  in  Figure 
6.1e,  which  represents  the  input  for  the  model.  The  corresponding  divergence  field  in  that 
figure  is  seen  to  retain  the  basic  north-south  convergence  line  in  the  proper  location,  but 
with  some  significant  fine-scale  differences  from  the  PAM-derived  divergence  field  in  Figure 
6.1b. 

The  divergence  field  computed  from  the  radar-derived  radial  winds  alone,  using  the 
method  of  Koecielny  et  al.  (1982),  is  shown  in  Figure  6. If.  Overlapping  sectors  for  this 
derivation  were  centered  on  every  other  grid  point  in  both  directions,  with  each  Cartesian 
“sector”  consisting  of  9  by  9  grid  points  representing  a  22.5  by  22.5  km  area.  This  single- 
doppler  derived  divergence  field  shows  general  agreement  with  those  obtained  from  the  PAM 
winds  (Figure  6.1b)  and  the  combined  radckr-radial  and  PAM-tangential  winds  (Figure  G.le), 
but  with  less  detail  and  with  lesser  extremes,  as  expected.  In  particular,  the  full  north- 
south  extent  of  the  convergence  line  is  diagnosed.  The  largest  differences  are  in  the  southeast 
portion  of  the  domain,  where  the  PAM  observations  are  sparse  and  the  Cartesian  “sectors” 
are  significantly  smaller  than  recommended  by  Koscielny  et  al.  (1982). 

B.  Case  Study  Simulations 

The  simulations  using  the  analyzed  Doppler  winds  were  initialized  with  special  soundings 
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taken  as  pairt  of  the  CINDE  program.  These  soundings  were  of  very  high  resolution  in  the 
vertical  (about  every  40-50  n?).  At  several  levels,  especially  in  the  boundary  layer,  the  winds 
were  of  questionable  reliability.  Therefore,  a  9  point  vertical  running  average  was  used  to 
filter  the  winds  before  the  models  were  initialized.  The  model  fields  were  interpolated  from 
these  soimdings  to  obtain  a  horizontally  homogeneous  initial  state  for  the  simulations. 

Starting  at  the  start  of  the  simulation,  the  low  level  model  winds  were  nudged  toward  the 
Doppler  winds  for  one  hour  with  a  nudging  timescale  of  one-half  hour.  After  one  hour,  the 
nudging  was  turned  off  to  simulate  an  actual  forecait  where  there  would  be  no  information 
available  on  the  wind  field.  The  mode'  was  then  run  from  one  to  two  hours. 

The  first  case  was  the  afternoon  of  17  July  1987.  A  north  to  south  low  level  convergence 
line  was  observed  to  develop  at  about  noon  and  remained  relaiively  stationary  until  convec¬ 
tion  developed  several  hours  later.  This  line  can  be  seen  in  the  PAM  and  radar  observations 
as  illustrated  in  Figure  6.1.  The  model  was  nudged  to  the  u  and  v  wind  components  of 
1315  LT  and  initialized  with  the  1430  LT  Denver  sounding.  The  vertical  motion  field  from 
the  full  model  version  at  1  hour  of  simulation  time  at  a  level  about  1  km  above  the  to¬ 
pography  is  shown  in  Figure  6.2a.  The  solid  contours  are  upward  vertical  motion  and  the 
dashed  contours  are  downward  vertical  motion.  A  region  of  upwcurd  motion,  showing  the 
convergence  line,  is  seen  Just  east  of  the  center  of  the  domain.  Consistent  with  the  generic 
nudging  simulations,  the  line  is  evident  in  the  nudged  model  fields  but  not  quite  as  well 
defined.  Compsoing  the  full  model  field  with  the  same  field  firom  the  mixed  layer  model 
(Figure  6.2b)  shows  that  the  main  upward  motion  areas  are  similarly  located  and  about  the 
same  magnitude  although  there  are  many  differences  in  the  fields. 

After  the  nudging  had  been  discontinued  for  1/2  hour  (at  1.5  hours  into  the  simulation), 
the  vertical  motion  fields  again  look  consistent  with  the  expectations  based  on  the  generic 
nudging  experiments.  The  upward  motion  is  replau:ed  by  subsidence  as  the  ‘*cold  dome” 
forced  by  the  thermally  imbalanced  upward  motion  collapsed.  Figures  6.3a  and  6.3b  show 
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very  similax  behaviors  for  the  full  model  and  the  mixed  layer  model,  respectively.  The  au-eas 
which  were  previously  rising  are  now  subsident  with  new  areas  of  rising  motion  forced  on 
the  edges  of  this  numerical  cold  pool.  By  the  end  of  the  simulation  (Figures  6.4a  and  6.4b), 
both  models  again  look  similar  as  the  fields  settle  down  after  the  shock  from  the  imbaJanced 
forcing. 

The  second  case  was  27  June  1987.  The  model  was  initialized  with  the  1300  LT  Denver 
sounding  amd  nudged  toward  the  radar/PAM  observations  at  1500  LT.  In  this  case,  con¬ 
vection  to  the  north  of  the  domain  produced  a  strong  outflow  boundauy  which  entered  the 
domain,  producing  an  east  to  west  convergence  line,  as  shown  by  the  PAM  observations 
in  Figure  6.Sa.  Figure  6.5b  shows  contours  of  the  east-west  wind  component  and  velocity 
vectors  at  1  km  above  the  ground  at  1  hour  of  simulation  time  from  the  full  model  version. 
The  corresponding  vertical  motion  field  is  shown  in  Figure  6.6a  where  a  region  of  upward 
motion  is  evident  at  the  convergence  line-  As  in  the  first  case,  the  mixed  layer  version 
was  very  similar  (Figure  6.6b).  After  the  nudging  has  been  turned  off  (Figures  6.7  and 
6.8),  again  similar  to  the  first  case,  the  upward  motion  decreases  by  1.5  hours  and  becomes 
subsident  by  2  hours. 

The  third  case  was  31  July  1987  and  was  substantially  more  complicated  than  the 
first  two.  On  this  day,  an  outfiow  boundary  from  convection  in  the  foothills  entered  the 
domain  and  collided  with  the  pseudo-Denver  convergence  line  and  other  gust  fronts  that 
were  already  in  the  domain  (Figure  6.9a).  Although  it  is  more  difficult  to  see  than  the  first 
two  cases,  the  models  behaved  similarly.  The  vertical  motion  fields  from  the  mixed  layer 
model  from  1,  1.5,  and  2  hours  are  shown  in  Figure  6.9b-d. 

To  siimmarize  these  real  data  cases,  the  behavior  of  the  models  in  these  cases  is  similar  to 
the  nudging  experiments  to  generic  flow  fields.  The  themudly  unbalanced  flow  forces  vertical 
motions  and  pressure  fields  which  attempt  to  oppose  the  flow.  As  soon  as  the  nudging  force 
is  discontinued,  as  would  have  to  be  done  for  a  real  time  forecast,  the  flow  quickly  responds 
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to  the  pressure  fields  and  becomes  very  different  from  the  observed  circulation. 

7  Summary  and  suggestions  for  future  research 

We  have  shown  in  preceding  sections  that  in  general  the  assimilation  of  just  wind  data  into 
a  dynamic  model  is  not  feasible  since  the  model  will  quickly  adjust  the  mass  field  in  such 
a  way  that  the  imposed  divergence  is  destroyed  by  perturbatioD  pressure  fields  that  oppose 
the  derived  divergence  fields.  The  situation  is  no  better  if  we  ha/f  a  system  of  measuring 
simply  the  msM  fields  on  the  mesoscale  as  the  model  will  develop  divergence  fields  that 
oppose  the  input  mass  fields.  What  is  clearly  needed  is  a  technique  for  observing  both  the 
wind  field  and  the  msiss  field  with  resolution  such  as  obtained  with  a  single  Doppler 'radar. 
Unfortunately  this  is  not  possible.  It  would  be  interesting  to  investigate  using  a  mix  of 
single  Doppler  radar  data  and  satellite  retrieval  of  the  mass  fields.  Current  satellites  can 
retrieve  soun<Ungs  with  a  horizontal  resdution  of  approximately  15km.  Vertical  resolution, 
however,  is  quite  crude.  It  is  possible  to  use  radar  data  to  enhance  the  vertical  resolution 
of  the  sounding  data.  The  radar  radial  velocity  data,  for  example,  can  define  the  top  of  a 
cold  pool  by  the  vertical  divergence  at  its  top.  By  assigning  the  depth  of  the  cold  pool  the 
satellite  retrieval  algorithms  can  be  modified  to  infer  the  average  temperature  of  the  cold 
pool.  Likewise  the  horizontal  boimdary  of  the  cold  pool  can  be  inferred  from  the  single 
Doppler  radar  data  yielding  a  more  precise  definition  of  the  boundary  between  the  cold  suid 
undisturbed  air  mass.  We  believe  that  this  is  a  natural  next  step  beyond  the  research  that 
we  have  accomplished  so  far.  Other  possibilities  of  obtaining  mass  and  wind  fields  using 
a  combination  of  remote  sensing  systems  may  also  exist  but  they  can  only  be  evaluated 
by  using  a  flexible,  modeling  system  such  as  RAMS.  Also  at  this  point  it  is  pre-mature 
to  restrict  four  dimensional  data  assimilation  multiple  data  systems  to  mini  computers. 
Once  a  satisfactory  technique  erf*  assimilating  a  mix  of  data  sources  in  a  mesoscale  model 
is  developed  using  what  ever  computer  resources  are  available,  then  the  technique  can  be 
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scaled  down  and  optimized  on  a  mini  computer  system.  VVe  clearly  have  not  reached  that 
level  of  maturity  of  data  assimilation  in  meso-beta-scale  models. 

Other  aspects  of  this  project  also  point  to  the  need  for  a  more  powerful  computer  than 
a  machine  rated  at  1  MIPS.  The  domain  size  of  100  km  by  100  km  is  far  too  small  for  a  6 
to  8  hour  forecast.  If  an  average  advective  speed  of  10  m/s  is  assumed,  the  interior  of  the 
domain  will  feel  the  effects  of  the  boundaries  within  an  hour  and  a  half.  The  forecast  period 
is  thus  limited  by  this  time  scale.  Forttmately,  computers  have  been  steadily  advancing  from 
a  price/performance  aspect,  so  now  a  machine  with  the  performance  about  20-40  times 
that  of  a  Micro  VAX  II  can  be  purchased  for  about  4  times  the  cost.  With  a  machine  of 
this  speed,  real-time  forecast  applications  of  the  full  model  version  can  be  realized.  For 
example,  a  coarse  grid  of  about  40  km  resolution  may  cover  an  area  the  size  of  a  state 
with  a  nested  grid  of  5-10  km  resolution  covering  a  local  area  of  interest.  The  effects  of  the 
model  initialization  may  then  be  removed  far  enough  from  the  region  of  interest  to  allow 
for  a  longer  forecast. 

The  utility  of  the  model  versions  developed  as  part  of  this  project  are  not  limited  to 
this  nudging  application.  Rather,  they  can  be  used  for  many  other  applications  where  their 
particular  features  are  desired. 
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figure  3.1.  Vertical  cross  section  of  the  potential  teirperature  field  used 
to  initialize  the  sea  breeze  test  siioulation. 


“igure  3.2a  Vertical  cross  section  of  the  potential  taopecature  field  in 
the  full  .nodal  version  of  the  sea  oreere  test  suisi  iation  4  hours  into  the 
sintilation. 
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Tifure  3.21)  vertical  cross  section  of  the  u-coeponent  of  n»tion  in 
tbi  full  Txsdel  version  of  the  sea  breeie  test  simulation  4  hours  into  the 
simulation . 


“igure  3.2c  Vertical  cross  section  of  the  vertical  motion  field  in 
cr.e  full  model  version  of  the  sea  oreeze  test  simulation  4  hours  in' 
simulation. 
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"igute  3. 2d  Vertical  cross  section  of  the  v-cotrponent  of  motion  in 
:he  fill  model  version  of  the  sea  oreeze  test  simulation  4  hours  into 
3  i.mxi^at  icr. . 
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"igur*  3.4b  Horizontal  cross  section  of  the  v-compor.ent  of  inotion  in 
t.~.e  ruxea-layer  rzsdei  version  of  t.-.e  .".eaceb  sloping  terrain  test  sionuiat 
1  ^-.ours  into  tbe  simulation. 
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Figure  3.4c  Horizontal  cross  section  of  the  w-component  of  motion  in 
the  mixed- layer  model  version  of  the  heated  sloping  terrain  test  simulation 
2  hours  into  the  simulation. 


iguE«  4.1a  Vertical  cross  section  of  the  u-con^nent  of  motion  in 
he  full  model  version  of  the  cartesian  wind  assimi*acion  test 
hour  i.'.to  the  simulation. 
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r-oure  4.10  vertical  cross  section  of  the  v-component  of  motion 
t.-.e  tu-1  .rnoael  version  of  t.-.e  cartesian  wind  assimilati.n  test 
1  hour  i.hto  the  simulation. 


tiqrur*  4 .  Id  Verticmi  cross  ssccion  of  ts*  p«rtucD«tion  pressure  field  in 
the  full  .Tiodel  version  of  tf'.e  cartesian  wind  assimilation  test 
I  hour  Into  the  si/taiiation. 
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Figure  4.1e  Vertical  cross  section  of  the  potential  ternperature  field  in 
the  full  model  version  of  the  cartesian  wind  assimilation  test 
1  hour  into  the  simulation. 
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"i<^re  4.2b  Vercicai  cross  section  o£  tbe  v-coitponenc  of  sotion  in 
tr.e  mixed-Layec  .rodei  version  of  tbe  cartesian  -ind  assimilation  test 
1  nour  into  tne  si-iajiacion , 
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Figure  4.2c  Vertical  cross  section  of  the  w-cowponent  of  nvotion  in 
the  mixed-layer  modal  version  of  the  cartesian  wind  assimilation  teat 
1  hour  into  the  aisulation. 
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Figure  4. 2d  Vertical  cross  section  of  the  potential  temperature  field 
t.he  mixed-layer  model  version  of  the  cartesian  wind  assimilation  teat 
1  hour  into  the  ainulation. 
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Figure  4.3a  Vemcai  cross  section  of  t.^e  u-cotnponenc  of  roti 
the  full  model  version  of  the  cartesian  wind  assimilation  test 
4  hours  into  tne  simulation. 
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Figure  <  .  4a  Vertical  cross  section  of  t>.e  u-component  of  motion  in 
tne  mixed-layer  model  version  of  the  cartesian  wind  assimilation  tes 
4  hours  into  the  simulation. 
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Figure  4.4b  Vertical  cross  section  of  the  w-component  of  motion  in 
-ne  mixed-layer  model  version  of  the  cartesian  wind  assimilation  test 
4  hours  into  the  simulation. 
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Figure  4.5c  Vertical  croaa  section  of  the  potential  teinperature  field 
the  full  model  version  of  the  cartesian  wind  assimilation  test 
5  hours  into  the  simulation. 
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isure  '  vemcai  cross  section  ot  the  u-coiteonent  of  siotion  in 
ne  full  T.odei  version  of  the  cartesian  wind  assimiiation  teat  for 
he  neutral  case  4  hours  into  the  siaaiiation. 
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“iqure  4.9  Vertical  cross  section  c£  the  u-conponent  cf  rroticn  in 
the  nujied-layer  model  version  ot  the  cartesian  wind  assimilation  test 
for  the  neutral  case  4  hours  into  the  simulaticn. 


igur«  4.^  Vertical  cross  section  of  the  u-corrponent  of  .totion 
he  full  .Tiodei  version  of  the  cartesian  ui.-.d  assimilation  test 
he  neutral  case  5  hours  into  the  simuiaticn. 
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Figur*  S.la  Horizontal  cross  section  of  t^e  u-conponent  of  niotion  in 
the  mixed'layer  model  version  of  the  radial  wind  assimilation  test 
1  hour  into  the  simulation. 
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^  '-"rizontal  tress  section  of  the  v-component  of  motion 

ed-layor  .mcoel  version  or  the  radial  wind  assimi.ation  test 
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Figure  5.1c  Horizontal  cross  section  of  the  w-component  of  motion  in 
the  mixed'*' layer  model  version  of  the  radial  wind  assimilation  test 
1  hour  into  the  simulation. 


rigur*  S.Ca  Horizontal  cross  section  of  ths  u-component  of  ration  in 
mixed-layer  model  version  of  tne  radial  wind  assimilation  test 
•)  hours  into  tne  simulation. 
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rigure  5.nb  Horizontal  cross  section  of  the  v-corrponent  of  motion  in 
ihe  mixed-layer  model  version  of  the  radial  wind  assimilation  test 
•»  hours  into  the  simulation. 
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Figure  5.2c  Horizontal  cross  section  of  the  w-con^ponent  of  motion  in 
the  mixed-layer  model  version  of  the  radial  wind  assimilation  test 
4  hours  into  the  simulation. 
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Figure  5.3a  Horizontal  cross  section  of  the  u-coeaponenr  of  motion  in 
the  mixed- layer  model  version  of  the  radial  wind  assimilation  teat 
6  hours  into  the  simulation. 
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Figure  5.3b  Horiicntal  cress 
the  mixed-layer  r.cGei  version 
6  hours  i.hto  t.he  simulation. 
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Figure  5.3c  Horizontal  cross  section  of  the  w-component  of  motion  in 
the  mixed-layer  model  version  of  the  radial  wind  assimilation  t 
6  hours  into  the  simulation. 
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Figure  6.1a  PAM  meaonet  observationa  of  wind,  temperature  and  dewpoint 
at  1500  local  time  (LT)  17  July  1987  (case  1) .  Bold  box  is  boundary  of 
100  by  100^.  km  computational  domain  (the  x  and  y  scales  differ),  large  dot 
denotes  the  off-centered  CP-3  radar  position  in  this  domain,  and  the  dashed 
line  A-B  denotes  the  dominant  north-south  convergence  line  on  this  afternoon. 
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Figure  6.1b  PAM  winds  from  Figure  6.1a,  objectively  analyzed  onto  the 
100  by  100  computational  grid  (with  an  extra  grid  point  in  each  direction 
to  accomodate  the  grid  staggering  used  in  the  model) .  Wind  vectors  are  shovm 
at  5  km  intervals  (every  other  grid  point) .  The  divergence  field  computed  from 
this  wind  field  is  contoured.  CP-3  radar  is  located  at  the  dot  near  the 
normalized  (x,y)  coordinate  of  (0.3, 0.5). 


NAOae-OCNIVCO  had l At  COMPOMCNT  IH/S) 


Figure  6.1c  Radar-derived  radial  velocity  field  over  the  computational 
domain.  Contoured  field  is  radial  wind  speed  (away  from  radar  is  positive) . 
Dot  denotes  the  radar  position. 
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Figure  6 . Id  Contoured  tangential  wind  component  of  the  PAM-derived  wind 
field  (cyclonic  motion  around  radar  is  positive) .  Dot  denotes  radar 
position. 


FtMt  U/V  0l«  iC-VSt 


X-NONN  1-33.9  TO  72.9  KHI 

Figure  6.1e  Final  vector  wind  field  for  mo<tel  input,  derived  by  combining 
the  radial  and  tangential  components  from  Figures  6.1c,d.  The  divergence 
field  computed  from  this  wind  field  is  contoured.  Dot  denotes  radar 
location. 
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Figure  6.1£  Divergence  field  as  derived  from  radial  wind  vector  field 
by  method  of  Koscielny  et  al.  (1982) .  Dot  denotes  radar  location. 
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r.rure  6.2a  Vertical  motion  contours  at  1  hour  simulation  time  from 
CINDE  case  1  at  1  'em  above  the  topography  from  t.he  full  model.  The  radar 
position  in  t.he  conputational  domain  in  this  and  all  subsequent  figures 
Is  at  (x,y)  coordinate  (-20. ,0.),  denoted  by  the  large  dot  in  this  figure. 
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figure  6.2b  Vertical  motion  contours  at  1  hour  simulation  time  from  Ci:«DE 
oase  1  at  1  <:n  above  the  topography  from  the  nuxed  layer  model. 
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Tiqure  6 . -lb  Same  as  **  ^  hours. 
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ri(jure  6.Sa  As  i.n  figure  6.1a  for  case  2,  at  ISIO  IT  2^  Jur.e  ISS"? . 
Outflow  boundary  from  shower  activity  north  of  t.he  domain  is  indicated. 
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qure  6.5b  J  toirronent  contours  and  velocity  vectors  at  . 

-Ti*  from  cr:rDE  case  2  ac  i  oi  aDov«  tJ'.®  copograpny  rrcm  tr.e  cu^a.  ••  i- 
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iqure  5.6a  Vertical  motion  contours  at  1  hour  stnniiatlon  time  from  CINDE 
ase  4  at  .  icn  above  tne  tcpograpny  from  the  full  model. 


I 


Vjam  MOOOi  t  X  Ok  &!•»  C*  1 


igure  i  .  i'o  Vertical  motion  contours  at  1  hour  simulation  time  from  C.NDE 
ase  at  1  <m  above  t.he  tcpograpny  from  the  muted  layer  model. 
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Figure  6.9a  .=i3  i.t  Figure  6.1a  for  case  3,  at  1430  LT  31  July  1987. 

Outflow  oounoary  from  mountain  a.oower  acti'/ity  west  of  the  cc.main  is  indicated 
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Figure  6.9b  '.'ectical  motion  contours  at  1  hour  simulation  ti-me  from 
OINTE  case  3  at  1  ion  aoove  the  topography  from  the  mixed  layer  model. 


Appendix  A 
THE  CSU  RAMS 


INTRODUCTION 

The  atimerical  atmospheric  models  developed  independently  under  the  direction  of 
William  R.  Cotton  and  Roger  A.  Pielke  have  recently  been  combined  into  the  CSU  Re¬ 
gional  Atmospheric  Modelling  System  (RAMS).  Development  of  many  of  the  physical  mod¬ 
ules  has  been  accomplished  over  the  past  15  years  and  has  involved  over  50  man  years 
of  effort.  RAMS  is  a  general  and  flexible  modelling  system  rather  than  a  single  purpose 
model.  For  example,  current  research  using  RAMS  includes  atmospheric  scales  ranging 
from  large  eddy  simulations  (A  s:  100  m)  to  mesoscale  simulations  of  convective  systems 
(Ax  as  lOOitem).  This  paper  will  discuss  the  options  available  in  RAMS,  the  engineering 
aspects  of  the  system  and  how  the  flexibility  is  attained. 

RAMS  OPTIONS 

RAMS  is  a  merging  of  basically  three  models  that  were  designed  to  simulate  different 
atmospheric  circulations.  These  were  a  nan-hydrostatic  cloud  model  (Tripoli  and  Cotton, 
1982)  and  two  hydrostatic  mesoscale  models  (Trimhack  et  aL,  1985  and  Mahrer  and  Pielke, 
1977).  The  capability  of  RAMS  was  recently  augmented  with  the  implementation  of  2- way 
interactive  grid  nesting.  Because  of  this,  the  modelling  system  contains  many  options  for 
various  physical  and  numerical  processes.  These  options  are  listed  below. 

The  following  options  are  currently  availstble  in  configuring  a  model: 

1.  Basic  ecpiations: 

Option  1  Non-hydrostatic  time-split  compressible  (Tripoli  and  Cotton,  1980) 
Option  2  Hydrostatic  incompressible  or  compressible  (Tiremback  et  al.,  1985) 

2.  Dimensionality:  1,  2,  or  3  spatial  dimensions 

3.  Vertical  coordinate: 

Option  1  Standard  cartesian 
Option  2  Sigm»-i 

4.  Horiiootal  coordinate: 

Option  1  Standard  cartesian 
Option  2  Latitude/longitude 

5.  Grid  Structure: 

•  Arakawa-C  grid  stagger 

•  Unlimited  number  of  nested  grids 

•  Unlimited  number  of  levels  of  nesting 

6.  Finite  differencing: 
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Option  1  leapfrog  on  long  timestep,  forward- backward  on  small  timestep,  2nd  or  4th 
order  flux  conservative  advection. 

Option  2  forward-backward  time  split,  6th  order  flux  conservative  (Turnback  et  ai, 
1987) 

7.  Turbulence  closure: 

Option  1  Smagorinsky-t3^e  eddy  viscosity  with  Hi  dependence 

Option  2  Level  2.5  type  clostire  using  eddy  viscosity  as  a  function  of  a  prognostic 
ttirbulent  kinetic  energy 

Option  3  O’Brien  profile  function  in  a  convective  boundary  layer  (Mahrer  and  Pielke, 
1977);  local  exchange  coefficient  in  a  stable  boundary  layer  (McNider,  1981). 

8.  Condensation 

Option  1  Grid  points  fully  saturated  or  unsattirated 

Option  2  No  condensation 

9.  Cloud  microphysics 

Option  1  Warm  rain  conversion  and  accretion  of  cloud  water  (r^)  to  raindrops  (r,), 
evaporation  and  sedimentation  (‘DipoU  and  Cotton,  1980) 

Option  2  Option  1  plus  specified  micleation  of  ice  crystals  (r{),  conversion  nucleation 
and  accretion  of  graupd  (r,),  growth  of  ice  crystals  (r^),  evaporation,  melting 
and  scdimentadon  (see  Cotton  et  oL,  1982) 

Option  3  Option  1  phis  option  2  plus  predicted  nucleation  and  sink  of  crystal  con¬ 
centration  (Ni),  conversion  and  growth  of  aggregates  (r^),  melting,  evaporation 
and  sedimenation  (see  accompanying  flow  diagram,  Figure  1).  The  nucleation 
model  includes:  sorption/depositicn,  contact  nucleation  by  Brownian  collision 
plus  thermophoresis  plus  diffiisiophoresis,  secondary  ice  crystal  production  by 
rime-splinter  mechanisnL 

Option  4  No  precipitation  processes 

10.  Radiation: 

Option  1  Shortwave  radiation  model  including  molecular  scattering,  absorption  of 
clear  air  (Yamamoto,  1962),  ozone  absorption  (Lads  and  Hansen,  1974)  and 
reflectance,  transmittance  and  absorptance  of  a  doud  la)rer  (Stephens,  1978), 
dear-cloudy  mixed  layer  approach  (Stephens,  1977) 

Option  2  Shortwave  radiation  modd  described  by  Mahrer  and  Pielke  (1977)  which 
indudes  the  effects  of  forward  Rayleigh  scattering  (Atwater  and  Brown,  1974), 
absorption  by  water  vapor  (McDosiald,  1960),  and  terrain  slope  (Kondrat’yev, 
1969). 

Option  3  Longwave  radiation  modd  induding  emissivity  of  a  dear  atmosphere  (Rodgers, 
1967),  emissivity  of  doud  layer  (Stephens,  1978).  and  emissivity  of  “dear  and 
doudy”  mixed  layer  (Herman  and  Goody,  1976) 
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Option  4  Longwave  radiation  model  described  by  Mahrer  and  Pielke  (1977)  includ¬ 
ing  emissivities  of  water  vapor  (Jacobs  et  oL,  1974)  and  carbon  dioxide  (Kon- 
drat’ycv,  1969)  and  the  computationally  efficient  technique  of  Sasamoh  (1972). 

Option  5  No  radiation 

11.  TVansport  and  diffiision  modules: 

Option  1  Advection-difiusion  model  (Segal  et  aL,  1980)  (To  be  implemented.) 

Option  2  Semi- stochastic  particle  model  for  point  and  line  sources  of  pollution  (Mc- 
Nider,  1981)  (To  be  implemented.) 

12.  Lower  botmdary: 

Option  1  Specified  surface  temperature  and  moisture  function  or  specified  surface 
fluxes  coupled  with  constant  flux  layer  condition  based  on  similarity  theory  (Man- 
ton  and  Cotton,  1977) 

Option  2  Surface  layer  temperature  and  moisture  fltixes  are  diagnosed  as  a  function 
of  the  ground  surface  temperature  derived  from  a  surface  energy  balance  (Mahrer 
and  Pielke,  1977).  The  energy  balance  includes  longwave  and  shortwave  radiative 
fluxes,  latent  and  sensible  heat  fliuces,  and  conduction  from  below  the  surfiKe. 
To  include  the  latter  effect,  a  multi-level  prognostic  soil  temperature  model  is 
computed. 

Option  3  Modified  form  of  Option  2  with  prognostic  surfrice  equatiosis  C&emback 
and  Kessler,  1985) 

Option  4  Same  as  Option  2,  except  vegetation  parameteruations  are  included  (Me- 
Cumber  and  PiellM,  1981;  1984;  McCumber,  1980)  (To  be  implemented) 

13.  Upper  boundary  conditions: 

Option  1  Rigid  lid 

Option  2  Rayleigh  Friction  layer  plus  Option  1-4 

Option  3  Prognostic  surface  pressure  (hydrostatic  only) 

Option  4  Material  surface  top.  (hydrostatic  osily)  (Mahrer  and  Pielke,  1977) 

Option  5  Gravity  wave  radiation  condition  (Klemp  and  Durran,  1983) 

14.  Lateral  boundary  oonditioos: 

Option  1  K1«wnp  and  Wilhelmson  (1978)  radiative  boundary  conditions 

Option  2  Orlanaki  (1976)  raffiative  boundary  conditions 

Option  3  Klemp  and  Lilly  (1978)  radiative  boundary  condition 

Option  4  Option  1,  2  or  3  coupled  with  Mesoscale  Compensation  Region  (MCR) 
described  by  Ikipoli  smd  Cotton  (1982)  with  fixed  conditions  at  MCR  boundary 
(see  Figure  2) 

Option  5  The  sponge  boundary  condition  of  Perkey  and  Kreitzberg  (1976)  when 
large  scale  data  is  available  frnm  objectively  analyzed  data  fields  or  a  larger  scale 
mc^el  nin.  This  condition  includes  a  viKous  region  and  the  introduction  of  the 
large  scale  fields  into  the  model  computations  near  the  lateral  botmdaries. 
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15.  Initialization 

Option  1  Horizontally  homogeneous. 

Option  2  Option  1  plus  variations  to  force  cloud  initiation. 

Option  3  NMC  data  and/or  soundings  objectively  analyzed  on  isentropic  surface 
wd  interpolated  to  the  model  grid. 

Option  4  NMC  data  interpolated  to  the  model  grid. 

As  one  can  see,  RAMS  is  quite  a  vex-satile  modelling  system.  RAMS  has  been  applied 
to  the  simulation  of  the  following  wesither  phenomena. 

1.  Towering  cumtili  and  their  modification 

2.  Mature  tropical  and  mid-latitude  cumulonimbi 

3.  Dry  mountain  slope  and  valley  drcul^ons 

4.  Orographic  cloud  formation 

5.  Marine  stratocumulus  clouds 

6.  Sea  breese  circulations 

7.  Mountain  wave  flow 

8.  Large  eddy  simxilation  of  power  plant  plume  dispersal 

9.  Urban  circulations 

10.  Lake  effect  storms 

11.  Tropical  and  mid-latitude  convective  systems 
ENGINEEaiNO  ASPECTS 

Because  of  the  large  number  of  options  in  RAMS,  the  structuring  of  the  code  needs  to 
be  careAiUy  considered.  This  section  will  discuss  various  aspects  of  the  code  structure  of 
the  system. 

Pre-processor  The  code  of  RAMS  is  written  in  as  close  to  the  FORTRAN  77  standard 
as  possible.  However,  with  a  program  as  large  as  this,  the  FORTRAN  standard  is  larking  in 
several  features  such  m  global  PARAMETER  and  COMMON  statements  and  conditional 
compilation.  To  remedy  these  insuflidencies,  the  RAMS  code  takes  advantage  of  a  pre¬ 
processor  written  as  part  of  the  RAMS  package.  This  pre-processor  itself  is  written  in  the 
77  standard  so  that  the  package  as  a  whole  is  highly  portable.  It  takes  full  advantage  of  the 
character  features  of  FORTRAN  and  has  executed  successfully  on  a  number  of  machines 
induding  a  VAX,  CRAY-1,  CRAY-X-MP,  and  CYBER  205  without  modification.  Some  of 
the  featmes  of  the  pre-processor  are  described  below; 
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1)  By  including  a  character  in  the  first  column  of  a  line  of  code,  that  line  can  be  “acti¬ 
vated”  or  “eliminated”  from  the  compile  file.  This  allows  for  conditional  compilation 
of  single  lines  or  entire  sections  of  code. 

2)  A  pre-processor  variable  can  be  set  to  a  value.  This  variable  can  then  be  used  in 
other  expressions  including  a  pre-processor  IF  or  block  IF  to  conditionally  set  other 
pre-processor  variables.  These  variables  also  can  be  converted  to  FORTRAN  PA¬ 
RAMETER  statements  which  can  be  inserted  anywhere  in  the  rest  of  the  code. 

3)  A  group  of  statements  can  be  delineated  as  a  “global”  which  then  can  be  inserted 
anywhere  in  the  code.  This  is  very  usefril  for  groups  of  COMMON  and  PARAMETER 
statements. 

4)  DO  loops  can  be  constructed  in  a  DO/ENDDO  syntax,  eliminating  the  need  for 
statement  labels  on  the  DO  loops. 

Two-way  interactive  grid  nesting  The  use  of  grid  nesting  allows  a  wider  range  of  mo¬ 
tion  scales  to  be  modeled  simultaneously  and  interactively.  It  can  greatly  ease  the  limi¬ 
tations  of  unnested  simulations  in  which  a  comprosnise  must  be  reached  between  covering 
an  adequately  large  spatial  domain  and  obtaining  sufficient  resolution  of  a  particular  locsJ 
phenomenon.  With  nesting,  RAMS  can  now  feasibly  model  mesoscale  circulations  in'a 
large  dnrm^n  where  low  resolution  is  adequate,  and  at  the  same  time  resolve  the  large  eddy 
structure  within  a  cumulus  cloud  in  a  subdomain  of  the  simulation. 

Nesting  in  RAMS  is  set  up  such  that  the  same  model  code  for  each  physical  process 
such  as  advection  is  used  for  each  grid.  This  makes  it  easy  for  any  desired  number  of  grids 
to  be  used  without  having  to  duplicate  code  for  each  one.  Also,  it  is  easy  to  add  or  remove 
a  nested  grid  in  time,  and  to  change  its  sixe  or  location.  There  is  still  the  flexibility  of 
choosing  many  modd  options  independently  for  different  grids. 

RAMS  has  adopted  the  two-way  interactive  nesting  procedure  described  in  Clark  and 
Fhrley  (1979).  This  algorithm  is  the  means  by  which  the  different  nested  grids  communicate 
with  each  other.  The  process  of  advancing  coarse  grid  A  and  One  nested  grid  B  forward  in 
time  one  step  begins  with  advancing  grid  A  alone  as  if  it  contained  no  nest  within.  The 
computed  Adds  from  A  are  than  interpolated  tri-cpiaciratically  to  the  boundary  points  of  B. 
The  inteHcr  of  B  is  then  cqrdated  under  the  influence  o!  its  interpolated  boundary  values. 
Finally,  the  fleld  values  of  A  in  the  region  where  B  exists  are  replaced  by  local  averages  from 
the  flelds  of  B.  An  increase  in  efficiency  over  the  Clark  and  Fhrley  method  was  implemented 
by  allowing  a  coarse  grid  to  be  run  at  a  longer  timestep  than  a  fine  grid. 

The  following  options  are  available  with  nesting  in  RAMS: 

1)  There  is  no  imposed  limit  (only  a  practical  one)  to  the  number  of  nested  grids  which 
can  be  used. 

2)  When  two  grids  B  and  C  are  nested  within  grid  A,  they  may  be  either  independent 
(occupying  different  space)  or  C  may  be  nested  within  B. 
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3)  The  increase  in  spatial  resolution  of  a  nested  grid  may  be  any  integer  multiple  of  its 
parent  grid  resolution.  Moreover,  this  multiple  may  be  specified  independently  for  the 
three  coordinate  directions. 

4)  A  nested  grid  may,  but  need  not,  start  from  the  groxmd  wd  extend  to  the  model 
domain  top. 

I/O  structure  For  those  machines  with  limited  central  memory  and  a  “non- virtual” 
operating  system  or  for  efficiency  on  virtual  systems,  RAMS  is  constructed  with  a  disk  I/O 
scheme.  When  the  scheme  is  operating,  a  subset  of  the  model’s  three-dimensional  variables 
will  reside  in  central  memory  at  any  one  time.  Computations  then  can  be  performed  with 
this  subset.  When  these  computations  are  finished,  a  new  subset  of  three-dimensional 
variables  are  requested  and  computations  performed  with  these.  The  RAMS  structure, 
thus,  is  dependent  on  this  I/O  scheme  and  consists  of  a  series  of  calls  to  the  I/O  scheme 
and  to  the  routines  which  do  the  calculations. 

Modiilarity  For  flexibility,  RAMS  is  written  as  modular  as  possible.  Eadi  individual 
physical  parameterisation  at  numerical  process  is  put  in  a  separate  subroutine  so  that  the 
routines  can  easily  be  replaced  for  different  options  or  with  new  developments. 

Computational  routines  The  routines  that  do  the  act\ial  computations  for  the  model 
are  written  so  the  implementor  of  a  new  or  replacement  routine  does  not  need  to 
be  concerned  with  most  of  the  details  of  the  rest  al  the  model  computations.  All  three- 
dimensional  variables  are  ’^passed"  to  the  subroutines  through  the  call  statement  with  other 
variables  passed  through  COMMON.  The  implementor  then  has  the  fledbility  to  structure 
his  routine  in  whatever  maimer  he  wishes  to  produce  the  desired  result.  This  concept  will 
uliirt  xnake  the  iziq>lexnentation  of  routines  from  other  models  and  programs  easier  with  less 
modification  required. 

Analyais  roun tines  A  set  of  subroutines  has  been  developed  for  analyzing  and  plotting  a 
variety  of  quantities  firom  fields  output  from  RAMS.  This  greatly  fadlitiates  the  interpreta¬ 
tion  and  imderstanding  of  modeled  atmospheric  phenomena.  The  quantities  diagnosed  by 
these  routines  include  vorticity,  divergence,  streamfunction,  energy,  momentum  flux,  most 
variances  and  covariances,  and  layer  averaged  quantities. 
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1  Introduction 

Four- dimensional  auimilation  of  obiervatioDal  data  into  a 
numerical  atmospheric  prediction  model  over  a  penod  prior 
to  the  start  of  the  pure  predictioo  (during  which  the  model 
receives  no  assistance  from  observatians)  it  ao  alternative  to 
a  single  3D  aatimilation  at  the  startup  time  of  the  pure  pre¬ 
diction.  In  iD  assimilation,  the  model  solution  is  weighted 
between  the  state  predicted  by  its  equations  and  the  state  de¬ 
fined  by  the  assimilated  data  so  that  some  informatiaQ  from 
earlier  data  inputs  is  preserved.  One  means  for  doing  this 
is  Co  ‘'nudge'  the  model  toward  the  input  data  by  mesms  of 
a  farcing  term  whose  strength  it  proportional  to  the  input 
data  minus  the  current  model  state.  The  success  of  AD  as- 
simUauon  depends  in  part  on  the  ability  of  the  model  fields 
to  properly  adjust  to  the  new  information.  The  assimilaced 
data  does  not  necessarily  represent  a  state  of  balance  either 
with  itself  or  with  the  current  model  solution.  This  is  par¬ 
ticularly  true  if  the  available  data  is  incomplete  in  the  sense 
that  it  contains  only  some  of  the  total  number  of  dynamic 
and  thermodynamic  fields  which  are  intimately  coupled.  For 
example,  there  are  currently  sridespread  plans  to  assimilate 
wind  information  from  2rE.XilAO  radars  and  sand  profilers 
into  operational  forecast  models  without  simultaneously  pro¬ 
viding  direct  observational  data  on  tonperature  or  pressure. 
It  IS  shown  in  the  present  study  that  this  can  have  a  very 
detrimental  effect  on  model  solutiosu.  Tsm  simulations  are 
run  m  which  the  h.jrixontal  winds  in  the  lowest  kilometer 
are  assumed  completely  known  in  real  time  by  means  of  re¬ 
mote  sensing  and  any  required  data  processing  ichems.  Ths 
honzontal  low  level  sands  in  the  numerical  model  a.';  midged 
toward  the  observed  sriod  field  to  attempt  to  bring  the  model 
state  into  closer  agreement  with  the  atmosphere.  Any  addi¬ 
tional  infonnation  oo  the  atmosphere  other  thsm  tbs  initial 
sounding  is  assumed  unavailable  to  the  model. 

2  Numerical  Experunents 

The  framework  for  combining  assimilated  sand  data  with  nu- 
mericai  integrations  of  the  governing  equations  is  the  Col¬ 
orado  State  University  Regional  Atmosphaic  .Modeling  Sys¬ 
tem  (RAMS).  It  is  run  in  its  ooohydrostatic  mode  with  20 
nonumfonnly  spaced  vertical  levels  extending  to  a  height  of 
I2fcm.  Ail  simulatians  are  conceived  es  having  slab  sym¬ 
metry  m  one  horisontal  d  'ction,  to  the  modal  n  run  as  a 


'fxgil  (4iD«  ft/fUiftiioo  <nch  Colofttdo  Slaiff  Loivwct 


2D  domAm  spaoning  I20krrx  in  the  other  horizontal  direction 
with  3jbm  resolution.  We  choose  the  following  simple  phys¬ 
ical  scenario  to  illustrate  the  main  point  of  this  paper.  A 
sounding  representing  a  honzontally  hotnogeneoua,  motion¬ 
less.  and  statically  stable  state  throughout  the  experimental 
domain  is  used  to  initialize  the  model.*  At  some  point  ax^er 
the  initialization,  a  cold  front  or  gust  &ont  moves  into  the  do¬ 
main  and  is  delected  by  Doppler  radar  which  is  fully  scanning 
the  lowest  kilometer  of  the  atmosphere.  The  front  is  observed 
as  a  6jtm  wide  convergence  line  advancing  at  4mj  *■  Into  tne 
still  air  and  being  trailed  by  Smj"^  winds  blowing  into  '.he 
front.  However,  no  mformatton  on  the  thermal  structure  oi 
the  front  is  available  to  the  model.  The  model  winds  are 
nudged  toward  the  observed  field  with  an  e-foiding  tune  of 
1800J. 

Figure  I  illustrates  the  perturbation  pressure  field  and  :he 
wind  vectors  in  the  plane  of  the  model  domain  five  minutes 
after  the  nudging  has  begun-  At  so  early  a  tune,  the  model 
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Fig.  1.  Vector  plot  of  the  2D  velocity  field  (peak  horror.: ai 
magnitude  u  l.Imj'M  and  contour  plot  of  the  perturbation 
presjure  field  (in  Paacal*)  for  the  stable  case  five  mmutes  a.-^er 
nudging  begins. 


winds  have  responded  directly  to  the  nudging,  all  other  forces 
.ncJuding  pressure  being  much  weaker,  and  a  oarrow  conver¬ 
gence  Line  begms  to  form  as  observed.  Following  a  persistent 
two-hour  application  of  the  nudging,  however,  the  structure 
of  the  model  winds  iihoim  in  fig.  2)  has  deviated  substa.-. 
•.aiiy  from  the  observed  structure  of  the  front  Lnsteao  .1 
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Fig.  2.  Same  as  fig.  1  but  two  hours  after  nudging  has  begun 
Peak  horizontal  velocity  is  S.-lrris"*. 


a  ccovergence  line  having  a  perpendicular  wind  component 
varying  Sms  ^  over  its  6km  width,  only  a  very  weak  and 
broad  zone  of  convergence  with  winds  varying  only  2ms"* 
a  100km  width  is  present .  The  reason  for  thU  is  that 
the  nudging  force  has  been  almost  completely  opposed  by  a 

horizontal  pressure  gradient  force,  which  is  apparent  in  hg. 
2.  The  surface  high  pressure  is  maintained  quasihydrostati- 
cally  in  the  vertical  by  a  dome  of  cold  air  (fig.  3)  formed  by 
upward  motion  of  stable  air  over  the  convergence  zone.  The 
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Fig.  3.  Plot  of  the  potential  temperature  field  coinciding 
with  the  fields  of  fig.  2. 


high  pressure  is  of  course  not  realistic  for  a  fronts  bound¬ 
ary.  In  summary,  the  modeled  acmosphenc  response  to  the 
nudging  of  the  velocity  field  was  to  counteract  it.  Had  the 
assimilated  data  also  contamed  the  infonnatioo  chat  the  air 
behind  the  front  was  significantly  colder  than  the  air  idiead  of 
it,  the  resulting  direct  nudgmg  of  the  model  temperature  field 
would  have  allowed  the  modeled  front  to  properly  develop  its 
convergence  line  structure. 


To  further  illustrate  the  problem  caused  by  tusimiiacing 
velocity  data  but  not  temperature  information,  we  cease  the 
velocity  assimilation  immediately  foliowing  the  two-hour  pe¬ 
riod,  as  if  the  pure  model  forecast  period  is  to  begm  then, 
and  continue  the  model  integration  with  no  further  nudg¬ 
mg.  Figure  4  shows  the  solution  one  half  hour  later.  V\Tsat 
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Fig.  4.  Same  as  fig.  2  but  one  half  hour  later,  nudging 
having  been  turned  off  during  that  half  hour.  Peak  horizontal 
veloucy  is  -S.Ottvi**. 


bad  been  a  convergence  zone  (and  still  should  be  according 
to  atmoapberic  observations)  is  now  a  region  of  divergence 
because  the  cold  dome  formerly  built  up  and  maintamed  by 
the  nudging  farces  is  in  the  process  of  collapsing.  This  solu¬ 
tion  is  obviously  grossly  out  of  agreement  srith  the  continuing 
frontal  passage  that  it  should  represent. 

The  underlying  fault  of  the  above  assimilation  experiment 
is  that  only  partial  obscrvaticnal  iafonnaiion  was  provided 
on  aa  entire  coupled  system.  In  the  example,  the  veloacy 
and  temperature  fields  are  coupled  because  (he  atmosphere 
is  stably  stratified.  In  a  nearly  neutral  atmosphere,  the  tem¬ 
perature  field  is  much  less  relevant  to  the  velocity  field,  and 
the  velocity  information  supplied  is  alone  tuAuent  to  pro¬ 
vide  a  good  solutio:^  nudging.  This  is  demonstrated  by 
a  second  expehmoit,  identical  to  the  first  except  that  the 
lowest  3km  of  the  atmosphere  are  very  cloae  to  neutral  in¬ 
stead  of  nearly  isothermal.  Figure  S  shows  the  veloaiy  fieid 
after  two  hours  of  nudging  to  have  a  sharp  convergence  line 
as  observed.  Since  the  upward  motion  and  compensating  di¬ 
vergence  above  ran  all  occur  within  the  deep  neutral  layer. 

no  significant  perturbations  of  the  temperature  field  are  gen¬ 
erated.  and  hence  no  large  fluctuations  in  the  pressure  heid 
result.  We  again  turn  off  the  nudging  as  in  the  first  exper¬ 
iment.  and  show  the  velocity  fields  one  half  hour  later  fig 
6).  .-Uthough  the  convergence  line  has  weakened  somewQat. 
It  IS  still  convergent.  Thus,  in  this  neutral  case,  nud^ng  has 
proved  a  benefit  to  the  model  solution. 
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Fi?.  5.  Same  ai  fig.  2  but  for  the  neutral  case.  Peak  hori¬ 
zontal  velocity  is  7.2mj'V 
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Fig.  6.  Same  a«  fig.  4  but  for  the  neutral  caee.  Peak  hon- 
lontal  velocity  is  6.4rTU"^ 


aloft,  the  formation  of  the  cold  dome  over  the  convergence 
would  have  been  realistic,  although  the  surface  high  would 
not.  One  might  also  consider  the  result  of  nudging  the  tem- 
peratxire  field  in  the  model  to  assimilated  data  through  a 
pseudo  heating  or  cooling  term,  without  any  velocity  obser¬ 
vations  being  available.  If  a  region  of  the  model  were  heated 
to  bring  It  closer  to  the  observed  temperature  of  a  warm 
atmosphenc  region,  and  that  atmospheric  region  had  been 
diabatically  heated  by  for  example  latent  heat  release,  radia¬ 
tion,  or  surface  heat  flux,  both  the  model  and  the  atmosphere 
would  tend  to  produce  upward  motion  in  the  heated  region 
and  agreement  would  be  good.  If,  on  the  other  hand,  -he 
observed  warm  region  of  the  atmosphere  had  instead  been 
heated  adiabatically  by  mechanically-forced  subsidence  tt.g. 
due  to  anticyclonic  vorticity  advection  aloft),  poor  agreement 
between  the  atmosphere  and  the  model  would  have  resulted. 

We  conclude  that  4Z5  assimilation  of  incomplete  attno- 
spheric  data  will  in  some  cases  improve  and  in  some  cases 
impair  the  solution  of  a  numerical  prediction  model.  .Vudg- 
ing  of  only  the  wmd  field  or  the  temperature  field  but  not 
both  to  observational  data  will  often  produce  an  unbalanced 
state  which  drives  the  model  solution  away  from  rather  than 
towju-d  the  observed  state.  Only  if  both  wind  velocity  and 
thermodynamic  information  can  be  assimilated  together  can 
AD  assimilation  be  expected  to  work  in  the  general  case. 
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3  Summary  and  Discussion 

This  study  has  illuscrated  an  example  in  which  <D  data  as- 
milation  into  a  numerical  prediction  model  produces  etro- 
.eouj  results  because  the  data  does  not  contain  complete, 
balanced  information.  In  this  case,  veloaty  information  is 
available  but  temperature  information  is  not.  FVoro  the  infor¬ 
mation  given,  the  only  rational  means  for  forcing  the  modeled 
velocity  field  to  the  observed  convergence  field  is  mechanically 
rather  chan  thermally.  Under  the  given  stable  atmosphenc 
conditions,  this  produces  an  anomalous  cold  dome  and  sur¬ 
face  high  pressure  at  the  convergence  Line  rather  than  the 
typical  fr'ontal  structure  assumed  to  accompany  the  conver¬ 
gence  line. 

It  is  interesting  to  consider  a  different  icenano  which  can 
hkewise  generate  surface  convergence.  If,  for  example,  ’he 
convergence  had  not  been  caused  by  baroclinity  at  the  sur¬ 
face.  but  instead  mechanically  by  cyclonic  vortiuty  advection 
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